

(NAS3L-CR-151960) DESIGN STODY OF A FEEDBACK 
CONTROL SYSTEM FOR THE MDLTICYCLIC FLAP 
SYSTEM FOTOR (NFS) Filial Report (Kamaa 
Aerospace Corp,) 90 p HC A05/MF A01 

CSCI 01A G3/02 


N77-24055 


Unclas 

29154 








,JI AEROSPACE 
Jim CORPORATION 

FIELD, CONNECTICUT 06002 


KAiUAK 

CORPORATION 

Report No. R-1494 
21 January 1977 


DESIGN STUDY OF A 
FEEDBACK CONTROL SYSTEM FOR THE 
MULTICYCLIC FLAP SYSTEM ROTOR (MFS) 

WASA rBISiSfcO 


FINAL REPORT 


1 Report No. 

NASA CR-151960 


2. Government Accession No. 


3. Recipient's Catalog No. 


4. Title and Subtitle 


pesign Study of a Feedback Control System for the 
Multi cyclic Flap System Rotor (MFS) 


5. Report Date 

21 January 1977 


6. Performing Organization Code 


7, Author(s) 

A. Weisbrich, R. Perley, H. Howes 


8. Performing Organization Report No. 

R-1494 


10. Work Unit No. 


9. Performing Organization Name and Address 

Kaman Aerospace Corporation 
91 d Windsor Road 
Bloomfield, Connecticut 06002 


11. Contract or Grant No. 

NAS2-8726 MODS 


12. Sponsoring Agency Name and Address 

1ASA Ames Research Center 
'•loffett Field, California 


13- Type of Report and Period Covered 

FINAL 


94035 


14. Sponsoring Agency Code 


15. Supplementary Notes 


'\na1ytTcal studies were performed to determine the feasibility of automatically 
providing higher harmonic control to a deflectable control flap at the tip of a 
helicopter rotor blade through feedback of selected independent parameters. 

Because the feedback study was based on the results of study work performed 
under Contract NAS2-7738, the rotor configuration data was updated to the con- 
figuration tested in the wind tunnel, which involved shortening the servoflap 
to half-size spanwise, and altering its balance and stiffness characteristics 

Following the updating of the analysis, control parameters were selected for 
input to the feedback system. The study work was then extended to arrive at a 
preliminary circuit design that would condition the selected parameters, weigh 
limiting factors and, ultimately, provide a proper output signal to the multi- 
cyclic control actuators. 

Since the first model of a feedback system would be tested on the wind tunnel 
multicyclic model, the preliminary design takes into consideration interfacing 
with existing multicyclic controls and uses existing test sensors. 


Study results indicate that feedback control for the higher harmonics is feas- 
ible; however, design for a flight system would require an extension of the 
present analysis. The present analysis was done for one flight condition - 120 
kts, 11,500 lbs gross weight and level flight - since a data base for other con- 




17. Key Words (Suggested by Author($)) 18. Distribution Statement 

Multicyclic Rotor Control 
Feedback System 
Controllable Twist Rotor 
Rotors, Helicopter 

19. S«cunty Ctassif. (of this report) 20. Security Classif. (of this page) 21. No. of Pages 22. Pt'C® 

UNCLASSIFIED UNCLASSIFIED 90 


*For sale by the National Technical Information Service, Springfield, Virginia 221 61 














KASAM 


A EROSPACE 
CORPORATION 


OLO WINDSOR ROAD. OLOOMFieLO. CONNECTICUT 06009 


lIAiUAS 

CORPORATION 

Report No. R-1494 
21 January 1977 


TABLE OF CONTENTS 


PAGE NO. 


LIST 

OF ILLUSTRATIONS 

iii 

LIST OF TABLES 

iv 

1.0 

INTRODUCTION 

1 

2.0 

MULTICYCLIC FLAP SYSTEM - ANALYSIS UPDATE 

2 


2.1 

Range of Study 

2 


2.2 

Analytical Procedure 

2 


2.3 

Multicyclic Flap Input 

5 


2.4 

Analysis Results and Discussion 

9 


2.5 

Concluding Remarks 

19 

3.0 

CONCEPTS OF MULTICYCLIC ADAPTIVE-CONTROL 

23 


3.1 

Control Function 

23 


3.2 

Optimization Parameter 

25 


3.3 

Rotor Characteristics 

27 


3.4 

Effects on Rotor Thrust Vector 

30 


3.5 

Search Strategies 

30 



3.5.1 Maximum Slope 

32 



3.5.2 Determination of Derivatives 

32 



3.5.3 Sequential Step Search 

33 



3.5.4 Magnitude of the Increment 

34 



3.5.5 Choice of Independent Variables 

35 

4.0 

SIMULATION STUDIES 

35 


4.1 

Effect of Various Optimization Functions 

36 


4.2 

Magnitude of the Increment 

38 


4.3 

Effect of Search Strategies 

40 


4.4 

Noise Investigation 

44 


4.5 

Dynamic Considerations 

47 


4.6 

Combined Noise and Dynamic Effects 

51 

5.0 

PRELIMINARY DESIGN 

51 


5.1 

Implementation Concepts 

51 


5.2 

System Test Sensors 

53 


5.3 

Processor Functions 

55 


5.4 

Modulator Function 

58 


i 




I[M\AK 

CORPORATION 

Report No. R-1494 
21 January 1977 


TABLE OF CONTENTS (continued) 

PAGE NO. 


6.0 CONCLUDING REMARKS 58 

6.1 Optimization Function Characteristics 58 

6.2 Characteristics of Measured Variables 58 ■ 

6.3 Control System Stability 60 

6.4 Rotor Characteristics 60 

REFERENCES 61 

APPENDICES 

A SURGEN PROGRAM A-1 

B BASIC PROGRAM, MCMODl B-1 

C MAXIMUM SLOPE C-1 

D SEQUENTIAL SEARCH MODEL D-1 


KAlWAX COI^PORATION 

OLD WINDSOR ROAD. BLOOMFIELD. CONNECTICUT OS003 


A 


OLD WINDSOR RDAD, BLOOMFieLO, CONNECTICUT OSOOR 



CORPORATION 


FIGURE NO. 

LIST OF ILLUSTRATIONS 

Report No. R-1494 
21 January 1977 

PAGE NO. 

la 

CTR Physical Properties 

3 

lb 

CTR Physical Properties 

4 

2 

CTR Control Optimization, ^ 

6 

3 

Optimization Functions 

26 

4 

Effect of Increment Size 

39 

5 

Response to Search Strategy 2D 

43 

6 

Effects of Iteration Interval 

50 

7 

Preliminary Design Concept 

54 

8 

Processor Control Functions 

56 

9 

Modulator Function 

59 

C-1 

Shortest Path Model 

C-2 

D-1 

Sequential Search Model 

D-2 


r 

i i .i 






A. 



CORPORATION 


Report No. R-1494 
21 January 1977 


LIST OF TABLES 

TABLE NO. PAGE NO. 

' 1 CTR Rotor Parameter Model Coefficients 10 

I 

2 Multicyclic Rotor Parameter Model- Coefficients 11 

3 Compilation of Multicyclic 6F Airloads Trim Cases 12 

4 SURGEN and 6F Trim Rotor Parameter Comparison (3315) 15 


5 SURGEN and 6F Trim Rotor Parameter, Comparison (3301X4) 16 

6 SURGEN and 6F Trim Rotor Parameter Comparison (330112) 17 

7 SURGEN and 6F Trim Rotor Parameter Comparison (3301X5) 18 


8 Effects of Higher Harmonic Control Input Phasing 

(PI and T3) 20 

9 Effects of Higher Harmonic Control Input Phasing 

. (Q7 and SI) 2.1 

10 Effects of Higher Harmonic Control Input Phasing 

(AW and L3) 22 

11 Surface Coefficients 28 

12 Derivative Coefficients 29 

13 Normalized Derivative Coefficients (Reordered) 31 

14 Effect of Various Optimization Functions 37 

15 Effects of Search Strategies 4T 

16 Effects of Noise 45 

17 Values of Derivatives 46 

18 Effects of Iteration Interval, t 49 

19 Combined Noise and Dynamic Effect 52 

C-1 Maximum Slope Model ' C-3 

D-1 Sequential Search Model D-3 


IV 




AEROSPACE 
CORPORATION 

OLO WINDSOR ROAD, BLOOMFIELD, CONHECTICUT 06002 

Report No. R-1494 
21 January 1977 



KMIAH 

CORPORATION 


1.0 INTRODUCTION 

Work performed under and in support of the Reference'’! Contract demonstrated 
that introducing a mixture of proper amplitude and phase of multiple harmonics 
(steady through 4Q) as a control command to a deflectable servoflap, located 
near the tip of a torsi onally-elastic rotor blade, could improve the airload 
distribution over the rotor disk. The results of this improved distribution 
are a delay in stall and compressibility effects, reduced hub vibratory shear 
loads, and reduced blade vibratory bending moments. Computer programs were 
developed by both the Contractor and by Ames Research Center that are effec- 
tive in predicting efficient control input settings for favorable tradeoffs in 
the important response parameters, given computer-calculated responses to an 
array of possible control combinations at a particular flight condition. At 
these efficient control settings, signifi.cant improvement over conventional 
rotor control was shown. 

Based on the results of the work discussed above, an existing experimental 
Controllable Twist Rotor (CTR), developed and wind tunnel-tested for USAAMRDL, 
Eustis Directorate, was modified to incorporate control modifications per- 
mitting the mixing of up to four harmonics of flap deflection, with indepen- 
dent control of the magnitude and phase of each. The control modifications 
and subsequent testing in the NASA Ames Research Center 40 x 80 wind tunnel 
were performed under Contract NAS 2-8726. Data from the wind tunnel tests are 
presently being analyzed for comparison to previous predictions. 

The next logical extension of the work performed on the manually-operated sys- 
tem, developed and tested under Contract NAS 2-8726, v/as to conduct a prelimi- 
nary design study to define and analyze a closed-loop self-optimization func- 
tion- The purpose of the work, which is the subject of this report, was to 
study the feasibility of automatically providing higher harmonic control 
through feedback of selected independent parameters. 

Because the feasibility study was to be based solely on the results of the 
initial study work performed under Contract NAS 2-7738, it was felt that the 
rotor configuration used in the study should be updated to the configuration 
tested which involved shortening the servoflap to half-size spanwise, and 
altering the mass balance and stiffness characteristics. The updated analysis 
is described herein. 

Following the updating of the analysis, control parameters were selected for 
input to the feedback system. The study work was then extended to arrive at a 
preliminary circuit design that would condition the selected parameters, 
weight limiting factors, and ultimately provide a proper output signal to 
the multi cyclic control actuators. Because the early study work was restricted 
to one flight condition (120 knot level flight), certain assumptions have been 
made concerning varying flight conditions and these assumptions are explained 
in the appropriate sections of this report. 


1 


COR^PORATION |[AiUAS 

01.0 WINDSOR ROAO. BLOOMFIELD. CONNECTICUT 06002 ' CORPORATION 

Report No. R-1494 
21 January 1977 


2.0 MULTICYCLIC FLAP SYSTEM (MFS) - ANALYSIS UPDATE 

2.1 Range of Study 

The primary purpose of the multi cyclic servoflap system is to reduce helicopter 
vibration levels by reducing rotor-generated vibratory loads that are trans- 
mitted to the fuselage. Although vibration problems are more severe at extreme 
operating conditions, the accurate prediction of rotor loads and performance 
at those conditions is difficult. Therefore, the prime thrust of this inves- 
tigation is in an area where retreating blade stall and high advancing blade 
Mach number effects are not significant. Rotor blade and disk loadings cor- 
respond to contempprary practice, and the propulsive force is representative 
of utility-type helicopters. 

The following flight condition and loading condition was investigated for the 
CTR with multi cyclic controls: 


Advance Ratio 

.333 

Disk Loading 

4.67 psf 

C^/o 

.092 

C^/a 

.0071 


The preceding parameters correspond to a rotor that has a diameter of 56 feet, 
a tip speed of 613 fps, and a solidity of 0.062. The ranges of disk loading, 
blade loading, and propulsive force loading correspond to the sea level flight 
conditions of a utility helicopter with a gross weight of 11,500 pounds and a 
flat plate drag area of 20 square feet. The study was conducted at an advance 
ratio of 0.333 (120 knots). 

Figures la and lb show the planform, inertia and stiffness of the rotor used 
in this study. 

2.2 Analytical Procedure 

The baseline rotor chosen for this study was the CTR developed and tested under 
the Reference 2 contract. This configuration differs from that examined in 
the Reference 1 contract in that the servoflap is reduced to half-size, and 
both mass balance and stiffness are increased. 

One disk loading was investigated with constant inflow across the disk 

(X = - .037). The rotor drag, (C /a was held constant at .0071, which is equi- 

valent to twenty square feet of flat plate area at y = .333 (120 knots). 

Lateral force, = 0. 

This rotor was optimized using the method outlined in Reference 3 for dual 
control rotor optimization. The following rotor parameters and limits were 
used as a measure of the effectiveness and determination of secondary control 
optimization: 
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a. Rotor horsepower ;< 750 

b. Maximum local blade angle of attack _< 12 degrees 

c. Maximum out-of-plane bending moment £ 9 in-kips 

d. Vibratory hub shears 5 220 lb. 

The angle of attack limit was chosen to provide a margin of maneuverability 
before local blade stall would be encountered. The 9000 in-lb peak-to-peak 
bending moment was selected on the basis of the calculated infinite blade life. 
The maximum out-of-plane moment occurs in the flap region of the blade where 
the endurance limit is +^4500 in-lb. 

Each of the four dependent variables outlined above can be expressed in quad- 
ratic form in terms of the independent control variables. These equations 
have linear coefficients and are developed for trim conditions only. Equa- 
tion (1) is the typical form of the relationship between a dependent variable, 
i.e., horsepower, and the collective and first-harmonic flap inputs: 

HP = ag + a^ + 82 6^5 + 83 6^^ + a^ 6 ^^ + a^ + ag 

+ 6q + ag 6^^ + 3g 6^^ 6^^ + a^g 

Similar equations were written for the remaining three parameters of interest. 
Appendix A is a listing of the regression analysis program SURGEN. The program 
uses a number of trim cases to calculate the rotor parameter equation coeffi- 
cients for the disk loading desired. In addition, the multiple correlation 
coefficients, standard error of the estimate, and a table of residuals are 
listed by SURGEN. Analysis of residuals for sum, sum of squares, mean, vari- 
ance and standard deviation, and a listing of the .orders of the residuals from 
the most negative to most positive are also tabulated for convenience. 

With the models, a tradeoff study was previously made to establish a region of 
flap control that would produce values of the four rotor parameters that meet 
the criteria established. Using a plot program on the Hewlett-Packard com- 
puter, representative contour plots, as in Figure 2, were generated. 

2.3 Multi cyclic Flap Input 

Multicyclic control was investigated with the same rotor for one gross weight 
corresponding to C^/a of .092. For C^/o of .092, the constant inflow across 

the disk was - .037. 

For the multicyclic flap concept, the range of cyclic control was limited to 
+ 5 degrees for each harmonic, such that the resultant maximum deflection for 
the 2nd harmonic and higher controls was 8 degrees. The flap ranges were 
selected based on the CTR work that had been conducted in the past. The 
restriction of the resultant of the 2nd and higher harmonic inputs to a maximum 
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120 kt 

y = .33 

11500 lb 

ZJa = .092 

907 lb 

ZJa = .007 

A 

613 fps 

6o = - 1 


Figure 2, CTR Control Optimization, 6 q = - 1. 
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of 8 degrees was to prevent excessive flap deflection from negating the bene- 
fits of higher harmonics. The following control range was established as a 
base: 


Steady 




«o = - 

V 

) 

First Harmonic 

Sine 

and 

Cosine 

6,3 f 3, 

, i 

he 

Second Harmonic 

Sine 

and Cosine 

+ 

5 

to 

Third Harmonic 

Sine 

and 

Cosine 

+ 

5 

to 

Fourth Harmonic 

Sine 

and Cosine 

+ 

5 

to 


Random Selection 


The selection of a control matrix by a method similar to that for the CTR was 
considered impractical because of the large number of control combinations pos- 
sible. For example: Involving only 2nd, 3rd and 4th harmonic controls, with 
seven control levels for each of six controls (i.e., 625* <^2c’ ^3s’ "^Sc’ ^4s’ 

a combined total of 117,649 cases are possible; A base of 140 cases 


with various combinations of 2nd, 3rd and 4th harmonics was selected by use- of 
a random number table. This method avoids individual prejudices and statisti- 
cally provides equal weight to the independent variables. Thus, a better 
chance is provided to account for the effects of each variable. within the lim- 
ited number of cases that were run. 


In Reference 3, cases trimmed at a fixed value of collective and first har- 
monic controls, together with variation in higher harmonic controls, showed 
horsepower, bending moment and maximum blade angle of attack to have signifi- 
cantly smaller changes from the base case values having no higher harmonic 
controls as compared to changes in vibratory hub shears (see Reference 3, 

Cases Al, 19, J3, M6, N8). Since the study objective was to arrive at optimum 
controls that reduce vibratory hub shears v/hile keeping bending moment, horse- 
power and maximum angle of attack with prescribed bounds, a single optimum set 
of collective and first harmonic controls was selected for the MFS study. 

Aside from the smaller variations of bending moment, horsepower and maximum 
angle of attack with higher harmonic control input, another reason for using a 
single set of collective and first harmonic flap control inputs is the very 
large number of trim cases which need to be generated to be able to create a 
reasonably accurate model which includes variations of collective and 1/rev 
controls. The number of linear regression model coefficients for each rotor 
parameter for the latter situation would have to be expanded from 27 to 54. 

The program scope did not permit such an expanded analysis. 


As previously stated, the set of collective and first harmonic controls chosen 
for case generation with higher harmonic control's was: 5^ = - 1 , = 3 and 

6-|^ = 5 (see shaded region of Figure 2i). Since variation of collective and. 

first harmonic input was eliminated, the linear regression mode] of each depen- 
dent rotor .parameter expands to: only 27 coefficients for the independent 
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control variables. Equation (2) is the form of the equation with hub shear 
being illustrated as one of the dependent variables: 



a-| ( 


"2 

'2c 

* ^3 hs 

U '3c 

^ "5 '4s 

^^6 

'4c 

+ a7 

2 

'2s 

+ 

Bg i 


2c 

• ag 

hs hs 

■ ®io '2s 

'3c ^11 

'2s 

'4s 



+ 

®12 


'4c 

+ a 

13 

^14 '2c 

'3s ^15 

'2c 

'3c 

"^16 

'2c '4s 

+ 

®17 


'4c 

+ a 

18 * 

®19 '3s 

'3c ®20 

'3s 

'4s 

+ 321 

'Ss '4c 

+ 

322 

^ 2 
^3c 

+ 5 

‘23 

he hs * 

®24 'sc 

'4c ^25 

'4s 

+ a 

26 '4s 

'4c 


^27 

« 2 
^4c 









(2) 

6 = 60 + 

'is 

sin 

^ + 

he 

cos if) + 6 

'2s 2 

4* *^2^ cos 2 

Ip + 

'3s ^ 

in 3 If) 

+ 

'sc 

cos 

3 ijj 

+ s 

sin 4 If) + 6^^ cos 4 if)’ 




(3) 


Due to the increased complexity of the higher harmonic flap inputs, it was not 
possible to use a graphic method to predict the best area of higher harmonic 
control inputs to minimize rotor parameters. Emphasis was placed on the use 
of successive models to predict trends. Since one purpose of this study was 
to determine the effect of higher harmonic control input on the reduction of 
4/rev pylon excitation, the major effort was oriented toward this goal. Since 
reduction of vibratory hub shears is of primary interest, the model for hub 
shears was used to predict higher harmonic controls which achieve minimum in- 
plane and out-of-plane vibratory hub shears. The model was introduced into 
the BASIC language program, MCMODl , Appendix B, which scanned all possible 
integer control combinations between - 4 degrees to + 4 degrees for each set 
of higher harmonic control (i.e., 62^, 62^} 632, 63^, 6^^, 5^^) and resulting 

in an absolute resultant higher harmonic control setting of less than or equal 
to 8 degrees. 


6F cases were then executed using selected model predicted controls that pre- 
dicted minimum out-of-plane hub shears. 
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2.4 Analysis Results and Discussion 

The following is a listing of the flight conditions examined for both the 
MCI Iti cyclic Controllable Twist Rotor (MCTR) and CTR: 

C^/a = .092 = 11500 lb rotor lift 

C^/o = .0071 Fy = 0 

p = .333 (120 kts) 

The CTR regression model (Equation 1) was based on about 42 trimmed cases. The 
model coefficients are tabulated for the rotor parameters of interest in Table 1. 

Figure 2, discussed previously and resulting from the above models, yielded 
good correlation between the models and actual trim. 

The MFS study was concerned with determining whether flap inputs with greater 
than first harmonic frequency could be used to reduce vibratory pylon loading 
without severe detriment to the remaining rotor parameters. 

Over 140 cases were trimmed at this rotor loading (11500 lb). This allows 
direct comparison with the CTR evaluated at this level. The fixed collective 
and first harmonic servoflap controls were chosen from the optimal CTR control 
region (6 q = - 1, 6^^ = 3, 6^^ = 5). 

The regression model coefficients generated by SURGEN for in-plane hub shears, 
out-of-plane hub shears, out-of-plane bending moments, horsepower and maximum 
angle of attack are presented in Table 2. 

Table 3 is a compilation of the multi cyclic trim cases generated showing 
values of the above rotor parameters for the control inputs presented. 

Table 4 presents the comparative values of in-plane hub shears, out-of-plane 
hub shears, bending moment, horsepower and maximum blade angle of attack for 
the multicyclic SURGEN model and the 6F trim program with only optimum collec- 
tive and first harmonic control input. 

These values are used for base comparison of cases having higher harmonic con- 
trol input. It is seen from Table 4 that, for the hub shears, bending moment 
and maximum angle of attack, the SURGEN model predicts high as compared to the 
6F program trim cases and, therefore, the SURGEN model appears conservative. 
Better agreement could be achieved by including more trim cases in the SURGEN 
model. However, for the purposes of this study, the model predictions are 
adequate. 

Tables 5 through 7 present cases having minimum out-of-plane hub shears within 
the limited trim cases generated. The out-of-plane hub shears for these latter 
cases are found to have reductions up to and over 70% with higher harmonic con- 
trol input. Bending moment and maximum angle of attack increased up to 27 % and 
horsepower values increased by as much as 7 % over the case values generated 
with no higher harmonic control input. 
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TABLE T. CTR ROTOR 

PARAMETER MODEL COEFFICIENTS 





HP 

BM . 

“max 

''v 


737.5 

7776 

13.38 

248.6 

"r 

0 

0 

0 

0 

^2 

2.367 

333 

.554 

- 11.68 

^3 

.777 

- 252 

- .564 

- 4.271 

®4 

0 

0 

0 

0 

^5 

.265 

28 

.028 

1.983 . 

"6 

.253 

18 

.066 

.453 

"7 

0 

0 

0 

0 


0 

0 

0 

0 

®9 

- . .123, 

13 

0 

- 1.125 

®10 

0 

0 

0 

0 

1 ZZ] 
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TABLE 2. MULTICYCLIC ROTOR PARAMETER MODEL COEFFICIENTS 



(Resultant) 

*^vop 

*^vip 


BM 


HP 


a 

max 

^0 

297.0 

287.32 

79.29 


4410.59 


744.94 


11.05 

^1 

- 17.59 

- 18.51 

- 5.83 

- 

142.22 


1.40 


.20 

®2 

9.34 

12.59 

- 6.13 

- 

206.07 

- 

1.25 


.59 


27.14 

25.17 

8.81 


32.57 

- 

2.92 

- 

.07 

®4 

1.62 

.07 

- .35 


52.89 

- 

1.25 


.01 


48.99 

48.86 

3.72 


10.43 

- 

.27 


.03 


- 15.76 

- 14,47 

- 13.51 


53.99 


.43 


.13 

®7 

2.17 

2.64 

.59 


37.94 


1.60 


.08 

^8 

8.84 

7.92 

1.45 


3.69 


1.12 

- 

.03 

ag 

2.49 

2.68 

.02 

- 

20.54 


.17 


.08 

®10 

- 2.77 

- 2.73 

.40 

- 

12.94 


.07 


.03 

®11 

.41 

1.07 

- .54 

- 

25.76 

- 

.57 


.03 

®12 

- 7.72 

- 7.92 

3.01 


22.42 


.18 

- 

.01 

®13 

12.58 

12,38 

2.70 


29.91 


1.50 


.12 

&14 

2.58 

2.79 

1.33 

- 

12.56 

- 

.06 

- 

.06 

®15 

1.30 

3.03 

- 6.29 

- 

49.34 

- 

.02 

- 

.09 

^16 

- 1.43 

.70 

- 4.48 

- 

5.04 

- 

.08 


.02 

^17 

■- 14.41 

- 13.82 

- 2.56 


35.07 

- 

.57 

- 

,06 

^18 

6.08 

6.39 

1.90 


60.62 


.78 


.04 

®19 

3.68 

4.25 

- .46 


16.38 


.25 

- 

.003 

®20 

16.39 

16.18 

1.16 


13.56 

- 

.03 

- 

,05 

®21 

- 6.94 

- 6.96 

- .94 

- 

98.58 

- 

.39 


.06 

^22 

11.06 

11.08 

3.79 


16.27 


.71 


.05 

^23 

13.11 

10.83 

8.50 


55,97 

- 

.05 

- 

.09 

®24 

24.10 

21.44 

9.09 

- 

12.51 


.36 

- 

.06 

®25 

15.13 

14.76 

3.10 


68.83 


1.36 


.07 

^26 

5.66 

6.12 

- 1.52 


30.84 

- 

.03 

- 

.006 

a?7 

19.60 

19.20 

3.15 


87.33 
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TABLE 3. COMPILATION OF HULTICYCLIC 6F AIRLOADS TRIM CASES (continued) 
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3.000 

— 5 ;ooo 

liOOO 

'“" 0 . 0 ' 
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0 . 0 
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0.0 

4.000 

0*0 

761. 0 

746 . 3 

146,6 

5700 . C 

764.0 

12 , 10 


oa 

-l.OOO 

3 * 3 ,/ 3 

5 . JJ^ 

0.0 

0.0 

3 . 536 

3.536 

0.707 

0*737 

787.0 

775.0 

160,9 

5807.0 

754,0 

U* 37 


09 

- 1.303 

3 . 3 J 0 

5,000 

0.0 

0.0 - 

2 . 121 

“ 2.121 

* 2.000 

0.0 

606*0 

6 DD .0 

115,4 

4391.0 

7 t 4.0 

n .69 


PI 

- 1.000 

3.000 

5.000 

1,414 

- 1.414 

1-414 

1.414 

0.3 

3*0 

346*0 

334.0 

108,2 

4913*0 

750,0 

11.41 


P 2 

- 1.000 

3.0 JO 

5.000 

0,707 

- 0 . 707 

1 . 000 

0.0 

4.000 

3.0 

dOL.o 

784*8 

159.0 

56 5 U 0 

766 . 0 

U .30 


— P 3 - 

'-moo ’ 

3,000 

- 5 iOOO 

'““ 4,000 

0.0 

0.0 

'' 0*0 " 

2.829 "* 

*“- 2 . 829 * ’■ 

768.0 * 

‘ 772*0 

158 . 6 

*'■ 5329*0 

" 7 G 5.0 

14 * 14 


P 4 

- 1.003 

3 . J 

5.000 

2 , 121 

2 . 12 L 

3.000 

0.0 

l. 4 t 4 

1*414 

678.0 

676.0 

116*8 ■ 

4183.0 

762.0 

12.62 


P 5 

-i;o 3 o 

3'.000 

5.000 

0.707 

- 0.707 

2.829 

2.829 

2.000 

0.0 

603.0 

778*4 

208,0 

5980-0 

756 . C 

11-16 


— P 6 *' 

- 1*000 

3 .J 33 

s.opo 

0-0 

0.0 

UOOO 

0-0 

2.121 

- 2.121 

708,0 

668*0 

170.2 

5147.0 

751 , 0 

11 . 80 


P 7 

-l.OOO 

3 . 0 UO 

5.000 

1.000 

0*0 

2.121 

- 2.121 

1.000 

0.0 

49 Z .0 

488.0 

112,0 

43 93.0 

7 -, 9,0 

11-24 


P 3 

- 1*333 

3 , JUJ 

5,003 

0.0 

0*0 

3.536 

3.536 

0.707 

3*737 

787.0 

775.0 

180.9 

5607.0 

754,0 

Uu 37 ■ 



- 1 , 000 — 

T.OJO' 

"'" 5.000 

0 * 70 T- 

“ 0.707 ‘ 

1.414 

1.414 ' 

0 * 707 “ 

— 0 ^ 737 ” 

” 511 .0 

* 506.3 

* 90,9 

4333,0 

' 746.0 

11.23 


Nl 

- 1*300 

3 , 00 j 

5.000 

uO ,707 

0.707 

1.414 

1 . 41 -* 

0.0 

0*0 

410.0 

403,4 

95,6 

44 J 3,0 

7 ‘* 4 , 0 

11 , 52 


N 2 

- 1.030 

3 ,‘JJO 

5,000 

2.121 

2*121 

0 . 707 

“ 0.707 

0.707 

“ 3.737 

554.0 

552.0 

99,0 

3707,0 

763.0 

18.56 


"'N 3 * 

- 1.300 

' 3.030 

5.000 

0-707 

0-707 

3.536 

‘ - 3.536 

2.000 

* 0.0 

‘ 757.0 ' 

752.0 

164.0 

4679.0 

762,0 

12.89 


NP 

-L .033 

3 .JJ 0 

5,000 

4,000 

0.0 

0.0 

0,0 

3.030 

3*0 

593.0 

580.7 

111.1 

4978.0 

775.0 

13 . 29 


N 6 

- 1,300 

3 .UJO 

5.000 

0*0 

0 . 0 

2 - 000 

0.0 

2.121 

- 2.121 

785.0 

764,0 

188,6 

5*6 i 0,0 

751.0 

11.55 


*' -Nr 

- 1.030 

3 * J JD 

" 5.000 

0 . 707 

0 . 707 

0 * 0 

0.0 


““ 3.0 ““ 

938.0 

' 92?- 6 

170 .? 

■ 5 >aoo. 0 

778 - Q 

13-1 J 


Nft 

-i*:o 3 

3 *uj 0 

5.000 

2.000 

0.0 

0.0 

0.0 

4.000 

3.0 

754.0 

7 J<'.d 

147*3 

5221,0 

7 lo. 0 

12.68 


01 

- 1,333 

3,030 

5.000 

2,121 

- 2.121 

0 * 707 

0.707 

l.ObO 

3 O 

372.0 

357.8 

105.0 

40 *, 6,0 

75 9.0 

11,43 


— 0 <;- 

- 1.300 

3.300 

5.000 

‘ 0,707 

‘- 0 . 7 O 7 

2.121 

“ 2.121 ' 

0.0 ‘ 

• 7.0 383.0 

• 376.7 

I 15 «l 

' 462 8,0 

750.0 

‘ 10*57 


07 

-l.OOO 

3.000 

5-000 

3.000 

0.0 

1. 414 

“ 1.414 

0-0 

3-3 

349.0 

348.9 

87.5 

4176.0 

7 b 3 * 0 

12 , 04 


08 

- 1.333 

3 . 3 jJ 

5 * 300 

0,707 

- 0 , 707 

3.536 

3.536 

0.0 

0.0 

647.0 

633*6 

192-5 

6217,0 

755*0 

11.85 


— 09 '" 

- 1*000 ' 

i.UOO 

5.000 

' 0.0 

- 0.0 * 

l.OOO 

‘ 0-0 

Z.lZl- 

” 2 ; 121 “* 

59&.0 

* 592*2 

67.0 

’ 5092,0 

765 , 0 

11.52 


SO 

- 1*000 

3.300 

5.000 

- 0.707 

0.707 

3*536 

3-536 

0.0 ", 

3.0 

669.0 

678,0 

172.0 

51 03 . 0 

751.0 

12 . 30 


SI 

- 1.330 

3.000 

5.000 

3,000 

• 0-0 

- 1.414 

1.414 

0.0 

3.0 

150.0 

147*0 

44,0 

4433,0 

768.0 

13*34 

ro :a 

ST’ 

- 1.000 V 

3 .UJO 

' 5.000 

- 0,707 

• 0 * 707 . 

“ 2.121 

0.0 ““ 

0.0 

3 .0 ““ 

29 2.0 

28 3.0 

70,0 

4722.0 

752.0 

U, 77 

*-» (D 

54 

- 1.030 

3 . JjO 

5,000 

- 2.121 

2.121 

0 . 707 

0.707 

1.000 

0.0 

569.0 

566 . 0 

116.0 

4773.0 

74 8.0 

12,57 

"O 

r*. 

S 5 

- 1*030 

3.000 

5.000 

0,0 

0.0 

2.000 

0.0 

“ 2.121 

2.121 

221.0 

208.0 

76,0 

5256,0 

76 0.0 

U. 9 B 

0) "S 

•“Ser- 

■-I.UOO “ 

3 , 0 J 3 T-' 

“S.OOO- 

'"O'; 707 

"" 0 . 707 “ 

- 3 . 536 “ 

3 V 5 J 6 *" 

“* 2 *. 700 — 0.0 

530.0 

“•* 504,0 

171.0 

5590,0 

• 772.0 

13,00 

rs r+ 

57 

- 1*300 

3 . 0 J 0 

5.000 

2,121 

2*121 

- 0. 707 

0.707 

“ 0*707 

0.707 

203,0 

' 193*0 

‘ 62-0 

4344,0 

771.0 

13 . 63 

C 

SB 

- 1*333 

3*000 

5.000 

1.000 

0*0 

- 2 . 121 

2.121 

I. 000 

3.0 

343.0 

322.0 

109.0 

4895,0 

759.0 

12.50 

Qi TZ 

■ — ■ — 



. . 

■ 

^ * ' — 



- — * 

■— 









— * TO 
CO I 


oi 

?*« 

jdS 

isfi> 


KAAIAK CORPORATION 

OLD WINDSOR ROAD, eLOOMFIELO, CONNECTICUT 06002 


TABLE 3 


COMPILATION OF MULTICYCLIC 6F AIRLOADS TRIM CASES (continued) 
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TABLE 4. SURGEN AND 6F TRIM ROTOR PARAMETER COMPARISON 





Case: 3315 Multicyclic CTR Without Higher Harmonic 

(MCTR w/o HH) 

Controls: = - 1 , = 3, 6^^ = 5 

Control Inputs 

ROTOR PARAMETERS 

SURGEN PREDICTION 

6F PROGRAM. TRIM VALUES 

(Resultant) 

297 

212 

D 

vop 

287.3 

210.3 

R . 

vip 

79.29 

44.8 

BM 

4410.6 

4175 

HP 

745 

748 

max 

11.05 

10.83 
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TABLE 5. SURGEN AND 6F TRIM ROTOR PARAMETER COMPARISON 


Case: 3301X4 

Controls: 5 q • 

Multi cyclic CTR, 
Out-of-Plane Hub 

= 1. 6^^ = 3, S2s 
•^Ic " ^2c 

Model -Predicted 
Shears 

= 4, 633 = 0. 6^ 
^3c ^ 

Controls for Minimum 
«4c = 1 

ROTOR 

SURGEN 

% VARIATION 
FROM 

6F PROGRAM 

% VARIATION 
FROM 

PARAMETER 

PREDICTION 

MCTR W/O HH 

TRIM VALUES 

MCTR W/O HH 

(Resultant) 

163.2 

- 45% 

118.7 

- 44% 

%op 

148.9 

- 48% 

113.7 

- 48% 

•^vip 

87.7 

10% 

36.3 

- 19% 

BM 

5344.6 

21% 

5340 

28% 

HP 

789.6 

6% 

791.7 

6% 

a 

max 

13.74 

24% 

13.78 

27% 
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TABLE 6. SURGEN AND 

6F TRIM ROTOR PARAMETER COMPARISON 



Case: 330112 

Multi cyclic CTR, 
Out-of-Plane Hub 

Model -Predicted 
Shears 

Controls for Minimum 

Controls; 6 q 

= - 1, <Si5 = 3, 6 

2s " *^3s " 

«4s = - 



6^,= 5, 6 

2c ^ ^3c "" ^ 

. «4c - 0 


ROTOR 

PARAMETER 

SURGEN 

PREDICTION 

% VARIATION 
FROM 

MCTR W/0 HH 

6F PROGRAM 
TRIM VALUES 

% VARIATION 
FROM 

MCTR W/0 HH 

(Resultant) 

143,6 

- 50% 

141.7 

- 33% 

^vop 

139 

- 52% 

140.0 

- 33% 

^vip 

68.7 

- 13% 

22.5 

- 50% 

BM 

5437.3 

21% 

5317 

1 

27% 

HP 

785.7 

5% 

786.2 

5% 

“max 

13.8 

25% 

13.82 

27% 

! 
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TABLE 7. SURGEN AND 6F TRIM ROTOR PARAMETER COMPARISON 


1 : • - : : " — i 

Case: 3301X5 

Multi cyclic CTR, 
Out-of-Plane Hub 

Model -Predicted 
Shears 

Controls for Minimum 

Control s : 6 q 

= - 1, = 3, 6 

2s " *^3s " " 

, 64s ■ 



= 5. 6 

2c " “ *^3c " ^ 

’ ^4c = 0 


ROTOR 

PARAMETER 

SURGEN 

PREDICTION 

% VARIATION 
FROM 

MCTR W/0 HH 

6F PROGRAM 
TRIM VALUES 

% VARIATION 
FROM 

MCTR W/0 HH 

(Resultant) 

156.7 

- 47% 

70.0 

- 67% 

*^vop 

148 

- 48% 

62 

- 71% 

*^vip 

75.8 

- 4% 

30.2 

- 33% 

BM 

5238.6 

19% 

5268.5 

26% 

HP 

781.8 

5% 

780.9 

4% 

%ax 

13.58 

23% 

13.58 

25% 
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In comparing the results of the SURGEN model predictions with the 6F trim cases, 
it is seen, as in Table 4, the SURGEN model is conservative. The SURGEN model 
predicts bending moments, horsepower and blade maximum angle of attack quite 
adequately. The hub shears are usually higher using the SURGEN model. However, 
the trends with higher harmonic control input are the same. 

It is also seen in Tables 5 through 7 that the controls that minimized the out- 
of-plane hub shears, in most cases, reduced the in-plane hub shears. When the 
in-plane hub shear SURGEN model was used to obtain minimum in-plane hub shears 
and these controls were then used in the out-of-plane hub shear SURGEN model, 
the out-of-plane hub shear increased 33%. Thus, it appears that the out-of- 
plane hub shear model is predominant, in that the out-of-plane hub shears are 
of the greatest magnitude and the controls that minimize the out-of-plane hub 
shears also reduce the in-plane hub shears. 

The SURGEN model’s limitations are possibly due to the relatively small number 
of trim data available (141 cases) upon which they are based. 

Variation in phasing of higher harmonic controls also shows pronounced effects 
on all rotor parameters. Tables 8, 9 and 10 present trim cases with variation 
in phasing of the 2nd, 3rd and 4th harmonic flap -control inputs, respectively, 
while amplitudes are held fixed. Limited phasing variations were used in the 
trim cases to generate the SURGEN models; therefore, optimum controls may not 
be obtained with 'the present models. 

Variation in optimum collective and first harmonic control input was not inves- 
tigated. 

2. 5 Concluding Remarks 

The following observations were made regarding the updated multicyclic analysis: 

• With the addition of 2/rev, 3/rev and 4/rev higher harmonic 
flap control input, reduction of the 4/ rev pylon excitation 
loads by at least 71% compared to the same rotor without 
higher harmonic control input was shown possible. 

Whereas the previous MFS study, cited in Reference 1, showed 
only 2/rev control input to be beneficial to reducing vibra- 
tory hub shears, the present study shows all higher harmonic 
controls together aid in reducing the vibratory hub shear. 

• Vibratory hub shears, bending moment, blade angle of attack 
and horsepower are significantly variant with both higher 
harmonic control amplitudes and higher harmonic control 
phasing. 
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TABLE 8. EFFECTS OF HIGHER HARMONIC CONTROL INPUT PHASING 



Case: PI and T3 

6F Trim, 2/Rev Phasing Effect 


Controls: = - 

1. = 3. 63^ = 1.414, . 0 



«lc = 5. 63, = 1.414, 6^^ = 0 



PI 

T3 

ROTOR PARAMETER 

60 = 2/135° 

62 = 2/-45° 

(Resultant) 

346 

456 

*^vop 

334 

450 

^vip 

108.2 

108 

BM 

4913 

4640 

HP 

750 

742 

a 

max 

11.41 

12.0 
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TABLE 9. EFFECTS OF HIGHER HARMONIC CONTROL INPUT PHASING 


Case: Q7 and SI 

Controls: “ “ 

6F Trim, 3/Rev Phasing Effect 
U <Sis = 3, ^22 = 3, 64^ = 0 
^Ic "" ^2c ^ ^4c "" ° 



Q7 

SI 

ROTOR PARAMETER 

63 = 2/135° 

0 

LO 

i 

II 

00 

*^v (Resultant) 

349 

150 

^vop 

348,9 

147 

^vip 

875 

44 

BM 

4175 

4433 

HP 

763 

768 

“max 

12.04 

13.34 
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TABLE 10. EFFECTS OF HIGHER HARMONIC CONTROL INPUT PHASING 


1 • ■ -H 

Case; AW and L3 

6 F Trim, 4/Rev Phasing Effect 


Controls: *Sg = - 

1, = 3, 625 = 4, 63 ^ = - 2 



^Ic " *^20 ^ ^ 



AW 

L3 

ROTOR PARAMETER 

64 = 2/-45° 

6 . = 2/135'^ 


^ 

(Resultant) 

354 

879 

^op 

344- 

864 

*^vip 

143 

180 

BM 

4616 

5845 • 

HP 

764 

762 

“max 

12.85 

12.03 
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9 Minimization of 4/rev out-of-plane vibratory hub shears show 
the most pronounced effect on overall reduction of vibratory 
hub shears. 

< f 

a In view of the pronounced effect of. phasing higher harmonic 
controls on rotor parameters j more comprehensive investigation 
of phasing is required, together with amplitude variation in 
higher harmonic" flap controls for trim case generation and 
subsequent rotor parameter SURGEN model development. 

• Incorporation of flap collective and first harmonic variation, 
at least within the optimum control region of the basic CTR, 
is anticipated to give a better understanding of the inter- 
relationship of all flap control inputs on rotor parameters 
of interest. 

3.0 CONCEPTS OF MULTICYCLIC ADAPTIVE-CONTROL 

3.1 Control Function 

The purpose of the multicyclic control is to minimi.ze certain parameters to the 
greatest extent possible while keeping others within acceptable limits. To 
accomplish this with feedback control, the definition of an "optimization" para- 
meter as a function of the various controlled parameters is required. Feedback 
is then used to vary certain controlling, or independent, variables in such a 
way as to minimize the optimization parameter. The controlled parameters that 
have been selected for this feasibility study are as follows: 

Bending Moment 

Horsepower 

Maximum Angle of Attack 
Hub Shear 

One of the criteria for selecting controlled parameters for the feasibility 
study was the ease with which they could be represented analytically on the 
basis of data from the 6F Aeroelastic Loads Program. When it comes time to 
implement this concept, a more important criterion will be the ease and repeat- 
ability with which the parameters can be measured on an operating rotor system. 
In this case, the control parameters are more likely to be: 

Control load related to bending moment 

Engine pressure 

Pitch link load 

Fuselage vibration levels 
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The controlling parameters, or independent variables, have been taken to be 
the sine and cosine components of the 2, 3 and- 4/ revolution servoflap control. 

The function of the feedback control system is to determine the effect of each 
independent variable (x's) oh the controlled parameters (y's) and, thus, on 
the optimization' parameter (P). The system then manipulates the x's to mini- 
mize P. Mathematically, the general case is: 

Optimization Function 


p = f(y-| 3 yg* y3> 

Change in P 


if ' I Ay-k 

k=l " 


where: 


^k ^2’ ^3’ ^4’ ^5’ ^6^ 


= I 


6 3y 


k 


4 

AP = I 
k=l 


i=l 

8P 


8X.. 


9yL 


6 3y-n 
I ^ Ax . 
i=l ‘^^i ^ 


The change in P caused by changes in a particular x^. is: 

AP - T SP , 

AP,- = I ^ • Ax- 

and the sensitivity of P to changes in a particular is: 


AP 


S. = 


1 _ 


3P 


3yi 


~N 


Ax. 


i k=l ^^k ^’^i 


From this, it can be seen that the feedback control system requires two kinds 
of information. The first is the relative importance of each of the controlled 
variables. The second is the sensitivity of each of the controjjled variables 
to changes in each of the independent variables. ^ The first is determined by 
the characteristics of the defined optimization function; and the second is 
determined by the rotor characteristics. 
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3. 2 Optimization Parameter 

The optimization parameter provides an integrated measure of the "goodness" of 
the rotor operation as described by measurements of .each of the controlled 
parameters. In effect, it provides a measure of the relative value of changes 
in each of the controlled parameters. An optimization parameter of the fol- 
lowing form has been used in the feasibility study; 


P = f^(BM) + f2{HP) + f3(AM) + f4(RV) 

The individual functions associated with each of the controlled parameters have 
been established on the basis of the operation of the feasibility study rotor 
with no multicyclic control and considerations of the consequences of each of 
these parameters in the operation of a helicopter. Figure 3 provides plots of 
the magnitudes and normalized slopes of each of these individual functions. 

The function for maximum angle of attack considers that the parameter is of no 
consequence as long as it is below the threshold value (12®). Operation above 
this threshold becomes undesirable at a very rapid rate. Bending moment and 
horsepower are similar to angle of attack, except that it is considered advan- 
tageous to reduce these parameters below their thresholds of 7000 in-lb and 
750 HP. In practice, it may be desirable to have the horsepower threshold be 
a function of some helicopter parameter such as forward speed; however, this 
feasibility study is only concerned with a single operating point of the rotor 
and, therefore, the threshold is fixed. Hub shear, and thus induced body vibra- 
tions, is considered to be of increasing concern from 0 on upwards. A decrease 
at a high level of hub shear is considered more valuable than the same decrease 
at a lower value of hub shear. The mathematical functions used are as follows; 


Bending Moment : 


BM < 7000 

f.| (BM) = 37 + 100 

BM > 7000 

f.|(BM) = 37 + 100 

Horsepower; 



HP < 750 

f2(HP) = 20 

HP 

750 " 

V, 

HP > 750 

f2(HP) = 20 

fHP 

750 " 


f BM 


7000 
f BM 


7000 


- 1 


- 1 


+ 10 


BM - 7000] 
1000 


- 1 


- 1+10 


HP - 7501^ 

^ J 


Maximum Angle of Attack : 


max 


< 12® 


> 12° 

- 10 (W - 12)“ 
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Hub Shear: 

2 ' 

f^iRv) = 100 . (§y 

Parameters such as horsepower and angle of attack could have been defined with 
discontinuities at the threshold points; however, operation of the feedback 
control system is better if the functions are continuous at all points. The 
scaling chosen in the functions above provides a value of the optimization 
parameter of about 100 at the rotor operating point corresponding to no multi - 
cyclic control. The measure of the improvement in rotor operation is provided 
by the decrease in this parameter from its initial value of 100. 

The use of this particular optimization function in this feasibility study is 
not a limitation on the applicability of the results. In fact, it was chosen 
because it represents a variety of subfunctions. The effect of some other 
optimization functions are discussed in paragraph 4.1. 


3.3 Rotor Characteristics 

The rotor used for the feasibility study was analyzed under a variety of multi- 
cyclic control conditions with the 6F Aeroelastic Loads Program, as described 
in Section 2 of this report. The data developed from this analysis were then 
used to fit a 7-dimensional quadratic surface for each of the controlled vari- 
ables. The coefficients for each of the four surfaces are shown on Table 11. 
The equation for each surface consists of a constant, six first order terms 
and twenty-one second order terms. The indices on the two dimensions of the 
matrix of coefficients represents the sine and cosine components of the second, 
third and fourth harmonic of the multicyclic control. For example, one term 
in the equation of the surface for bending moment is: 


- 12.94 . X23 


^3c 


The dimension used for the x's is degrees. 


The real interest is in the slopes of these surfaces as represented by the 
derivatives of each of the controlled parameters with respect to each of the 
independent parameters. The coefficients of the derivative equations are given 
on Table 12. The equation for the derivative of a controlled variable with 
respect to an independent variable contains one first order term and six second 
order terms. The coefficients for the first order term and the off-diagonal 
second order terms are the same as the coefficients for the surfaces; however, 
the coefficients for the second order terms on the diagonal are double the value 

2 

for the surface since they correspond to an x term. 

An inspection of these coefficients gives little, if any, indication of one 
independent variable being dominant with respect to any controlled variable. 

For instance, while the second harmonic would dominate bending moment when 
there is little or no multicyclic control (all second order terms at or 'near 
0), 3s, 4s or 4c components with magnitudes in excess of 2° would have a 
greater effect, as would 3s in combination with 2° of 4c, or vice versa. 
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TABLE 11. SURFACE COEFFICIENTS 



• FIRST.' !v 

SECOND 

ORDER 

' 



t- 

ORDER;, 2S 

2C • 

3S- 

3C 

4S 

AC : 

COEFFICIENTS FOR BMJ i 

CONSTANT^ 

4410. S9i 



2S 

-;l 42.220 37.94X) 

3.690 

”20.540 

-12.940 • 

"25.760 

-22*42,0 

2C> 

••••206 .070 

29.910 

-12.560 

-49.340 

"5.040 

35 ♦ 070 

33' 

32 1 570 


60.620 

16.380 ' 

"13.560 

-98.580 

3C 

52.890 , ;• 



16.27'’0 

55.970 

-12.510 

4S" 

-10.430 


/ 


68 ♦ 830 

30.940 

4C ' 

53 .'990 





87.330 

COEFFICIENTS FOR HPJ 

CONSTANT- 

744.94 



- 

2S- 

l'.4p0 I.600 

1.120 

0.170 

0.071 

-0.570 

0.180 

2C‘ 

-1-.250 

1.500 

-0.059 

-0.019 

-0.082 

-0.570 

3 S.- ■ 

-2.920 

^ y 

. 0.780 

0.250; 

-0 ♦ 033 

-0^90 

ZQ'-' 

-1-.250 



0.710 

-0.050 

0.360 

43 

•^a .270 , > 




1.360 

-0.030 

4C 

0.430 

• 


. 


1.190 

COEFFICIENTS FOR AM-J 

CONSTANT-:^* 

1 1.053 



■ ; 

4 m Wj 

0.V196- 0.O79 

-0.027‘ 

“0 . 082 

0.030 

•^0.030 

-0.015 

2C- 

• 0.594 . • ‘ 

. f0*125 

-0.061 

. -0.087 

0.619 

.-0i057 

3S. 

. -OfOiSO ' 


0.041 

; -0.003 

-0 . 054 

0.063 

3C- 

-0..0'09 Y 

«• » 


■ 0.04,7 

-0 . 006 

-0.059 

4S 

-a. 032 ' 


* 

■/ 

0.071 

-0.006; 


“0.130 •. 

. * 

' * r 

*/ \ 
r. 


f , 

•’ 0 . 059 ; 

COEFFICIENTS FOR RV4' 

CONSTANT* 

2'97-'/' 

i 


* s * 

2S. 

“i7;^5.9o 2V170 

8 ♦'840 

. 2.4930 

'i- 2 . 77 ‘ 0 - 

’0.410 

-y\720 

2C 

9.340 

12. *580,, 

2.5,80 

1.'300 

-1 . 430 

“14.V410. 

3S‘; 

27. ,140 


6 

,'• 3.680 

16.590 

“6 -.940 

3C ‘ • 

■ lv'620 


S'. 

■ ll.;060“ 

13.110 

24.io'o; 

43 ' 

•*48;.990 




15.130 

5 ♦ 660’ 

4C.' 

■ -15.t760 



' 

t 

19,.6'00 

•a * 

DONE 

t 
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TABLE 12. DERIVATIVE COEFFICIENTS 


FIRST SECQNIi ORDER 



ORDER 

2S 

or 

Am K / 

3S 

3C ' 

4S 

4C 

COEFF'ICIENTS FOR DERIVATIVES OF 

BMJ . 


. 


2S 

"142 *220 

75 * 880 

3*690 - 

■20*540 

--12*940 

•~25 * ’76'0 

"22*420 

2C 

"206 ♦070 

3*690 

59*820 •• 

12.560 

"49,340 

•"5 * 040 

35 * 070 

3S 

32*570 

••■20*540 

"12*560 :i 

L21.240 

16*380 

"13*560 

•• 98. 580 

30 

52*890 

• J 2*940 

"49*340 

16*380 

. 32*540 

55 <970 

"12.5;i0 

AS ' 

"10*430 

"25*760 

"5*040 • 

-13*560 

55*9/0 

137*660 

30*940 

AC 

53*990 

"22 * 420 

35.070 - 

-98*580 

-12*510 

30*940 

174*660 

COEFFICIENTS FOR DERIVATIVES OF 

HP *' 




2S 

1*400 

3.200 

1.120 . 

0.170 

0.071 

"0*5/0 

0*180 

2C 

"1*250 

1*120 

3*000 

"0*059- 

■ "0*019 

"0*082 

"0*570 

3S 

' "2*920 

0*170 

"0*059 

1*560 

0^250 

"0 * 033 

"0>390 

3C 

"1 * 250 

• 0*071 

"•0*019 

0*250 

' 1 * 420 

"0*050 

0 * 360 

4S . 

"0 * 270 

•0*570 

-0*082 

"0*033 

0*050 

2*720 

"0 * 030 

4C- 

0*4i^0 

0*180 

"0*570 

-0*390 

0*360 

"0 * 030 

2 <300 

COEFF 

ICIENT-S FOR DERIVATIVES OF 

MM 




2S 

0*196 

0vl58 

"0*027 

-0*082 

0 ♦ 030 

"0 * 030 

-0*015 

2C 

0*594 

••••0*027 

, 0*250 

"0*061 

-0i087 

0*019 

-0 * 057 

3S 

"0*068 

"0*082 

"0.061- 

0 *'082 

-0 * 003 

-0 * 054 

0 * 063 

3C 

"0*009 

0i030 

r0*0S7 

"0*003 

0*094 

-0*086 

"0*059 

AB- 

"0*'032 

"0*030 

0*019. 

"0*054 

■ "0*086 

0*142 

-0*006 

AC 

-O’* 130 

"0*015 

"0*057 

0*063 

-0*059 

"0*006 

0*117 

COEFFICIENTS F 

OR DERIVATIVES OF 

RVJ 

» 



2Sv 

••••17*590 

4 .-340 

. 8^840 

2 *<490 

-2*770 

0*410 

"7*’;^20 

2C 

9*340 

8*840 

25*. 160 

2*580 

, , 1*300- 

"1*430 

- 14*410 

3S' 

27*140 

2*490 

■= 2*580 

12*160 

'• 3*680 

16*390 

"6*940 

3C • 

1 * 620 

"2*770. 

1 * 300 ‘ 

3 * 680 

’22*120 

13* no 

24*100 

4S 

48*?^.90 

0*410 

"1*430‘ 

16*390 

13*110 

30*260 

5*660 

AC 

"15 '♦•760 

"7*720 

-:14*410! 

"6*940 

\ 24 #100 

5*660 

39*200 


HONE 
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Another view of these coefficients is provided on Table 13. In this case, the 
.coefficients have been reordered to bring together all of the derivatives with 
respect to a particular independent variable. In addition, the coefficients 
have been normalized to 1% of the threshold value used in the optimization 
function for the corresponding controlled variable. The arrangement of this 
table provides a view of the effect of each independent variable on the opti- 
mization parameter; however, for the picture to be complete, each of the coef- 
ficients should be multiplied by the slope of the corresponding optimization 
function for a particular operating point as given in Figure 3. This indicates 
that each of the independent variables has significant impact on the optimi- 
zation function, even the 3c component which has small first order coeffi- 
cients, has second order coefficients associated with 3c, 4s and 4c which are 
sizeable. 

While it is dangerous to make sweeping conclusions based on rotor performance 
data at a single operating point, there are two conclusions that can be drawn 
that are important to the remainder of this feasibility study. First, the 
second order coefficients of the derivatives can be dominant under a variety 
of operating conditions. This tends to confirm the necessity for determining 
the derivatives during actual operation of the rotor, rather than incorporating 
them as part of the pre-programmed control function,. The second conclusion is 
that each of the six independent variables has significant effect on at least 
bending moment and shear and, therefore, must be incorporated as part of the 
feedback control function. 

3.4 Effects on Rotor Thrust Vector 

In general, the changes -in multicyclic control introduced at the tip end of 
the blade to minimize the controlled parameters also change the lift and hori- 
zontal thrust vector developed by the rotor. Therefore, to maintain rotor 
trim, it would be necessary to introduce changes to the 1/rev cyclic pitch at 
the root end of the rotor. If this were not done and if the effects of the 
multicyclic control were large, the changes to multicyclic control would cause 
changes in helicopter operation which would, in turn, require a constant pilot 
attention in order to maintain the required operating condition. However, 
work with the 6F Aeroelastic Loads Program has indicated that the effects of 
multicyclic control on rotor trim are, in fact, small at the operating point 
investigated. Therefore, for the purposes of this feasibility study, the 
effects of multicyclic control on rotor trim are taken as not large in compari- 
son to other effects for which the pilot must continuously compensate. This 
assumption would probably not be valid if the rotor tip end pitch changes at 
1/rev were made subject to feedback control. 

3.5 Search Strategies 

The feedback control system requires a strategy to direct its search to find 
the path from the current operating point to the "best" operating point as 
defined by the optimization parameter. This strategy is contained in the 
logic that is used to analyze the measurements from the operating rotor in con- 
junction with the characteristics of the optimization function to determine 
which of the independent variables to increment next. 
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TABLE 13. NORMALIZED DERIVATIVE COEFFICIENTS (REORDERED) 



I-.LRST 

SECOMB 

ORDER 





ORDER 

2S 

2l‘ 

3S 

3C 

4S 

4(.' 

i.tQ\M;oAr’;i:uES‘ wm-i 

RESPECT 

TO 2S 





0M 

■~2<0Z2 

t,084 

0,053 

"0.293 

”0, 105 

"0<368 

"0 . 320 

iip 

0,18/ 

0,42’/ 

0,149 

0,023 

0,009 

"0,076 

0.024 

AM 

1 , 635 

1,J20 

-0,-228 

•0,633 

0,252 

"0 , 250 

••••0,124 


•“5,063 

1,447 

2,947 

0 . 830 

•■•0,923 

0 ^ 1 37 

-2,5/3 

DERigAT.].OEG WITM 

RESPECT 

TO ,2C 





BM' 

•■•2,944 

0 .. 053 

0,855 

“0,179 

“0, /05 

“0,072 

0,501 

HP 

•■•0,167 

0,1.49 

0 , 4p0 

- 0,008 

•0,003 

•••0,0 'Ll 

“0,076 

AM 

4 , 953 

0,220 

2 , 082 

"0.51 1 

•••o;7?6 

I 0,158 

• 0.474 

RP 

3,113 - 

2, 9-r/ 

8 , 387 

0,060 

0 , 433 

•0,477 

".4 , 803 

riERlVATJVES w.mi 

R(;:spr.(.n 

rd 3S 





Bit 

0 , 465 

■0>293 • 

-0,179 

.1 , 732 

0 , 234 

"0,194 

'" I. ,408 

Hl=- 

-0,'309 

0,023 

•0,008 

0v208 

0,033 

"0,004 

-•0,052 

AM 

-0,565 

0,603 

■••0,5.1 .1 

0,687 

■0,026 

"0.447 

0,524 

RkJ 

9,047 

0,330 

0 , 860 

4,053 

1,227 

5, 463 

“2,313 

riERiVATiOEs uu;th 

RESPECT 

fO 3C 





BM 

0 , 756 

•0. 105 

■0,705 

0,234 

0,465 

0,000 

•”0-,179 

IIP 

“0,167 

0,009 

“0 , 003 

0 , 033 

0,189 

“O,, 007 

0,048 

AM 

•“0 , 077 

0 , 252 

•0 , 726 

"0,026 

• 0,707 

"6,718 

"0,495 

RU 

0,540 

0,923 

0,433 

1,227 

7,373 

4,370 

8,033 

DERI VAT ;i.VES WCTH 

RESPECT 

rn 48 

V 




BM 

“0,149 

0,?6S 

“0<(y72 

"0,194 

0 ♦ 800 

1,967 

0,442 

HP 

■••••0,036 

0,076 

"0 , 01 . 1 . 

"0,004 

"0.007 

0 , 363 

"0,004 

AM 

“0,263 

-0,250 

0,158' 

"0,447 

"0,718 

1,180 

"0 , 05.1 

RV 

16, '330 

0,;I37 

”0,477 

5 ,463 

4,370 

10,087 

1,887 

DERIVATIVES- WITH 

RESPECT 

TO 4'C 





BM 

0;’77J. 

•0 ♦ 320 

0,501 

•“.I. ,408 

"0‘,179 

0,442- 

2, '4 95 

HP 

0,057 

0,024 

••••0,076 

"0 , 052 

' 0,048 

"0.004 

0,317 

AM 

"1,085 

•0,124 

•0,474 

0,524 

"0 , 495 

"0,051 

0,978 

RV 


-2, "573 

“ 4 , 803 

"2,313 

8 , 0.^3 

1 , 88-7 

13,067 
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3.5.1 Maximum Slope . Theoretically, the shortest path from a given operating 
point to the "best" operating point is found by moving at all times in the 
direction of maximum slope on the seven-dimensional surface relating the six 
independent variables to the optimization parameter.. Movement continues until 
the slopes in all directions become. zero. In a practical system where move- 
ment takes place with a finite increment, the point where all slopes equal 
zero is not found. Instead, the slope changes sign as. the minimum point is 
passed. The ability to. get closer to the minimum point, .is inversely propor- 
tional to the size of the increment used. On the other hand, the time required 
to move from one operating point to another is also inversely proportional to 
the size of the increment. Therefore, a tradeoff or provision of a variable 
increment. size is required. However, in a practical system, the noise associ- 
ated with the measurements of the operation of the rotor also limits the 
degree to which the minimum point can be approached. 

To implement this search strategy, the control system must first determine the 
derivative of the optimization parameter with respect to each of the indepen- 
dent variables. These derivatives are then inspected to select the ;one with 
the maximum magnitude. The independent variable corresponding to that deriva- 
tive is then incremented in a- direction corresponding to. .the sign of the deriv- 
ative. 




3.5.2 Determination of Derivatives . As indicated in paragraph 3.1, the con- 
trol system must be responsive to the derivative of the optimization function 
with respect to each of the independent variables. This may be expressed as 
the sum of the products .of the derivative of the optimization function with 
respect to the controlled variable and the derivative of the controlled vari- 
ables with respect to the independent variable. Knowledge of the first deriv- 
ative comes from the optimization function specified. Knowledge of the other 
derivatives must come from measurements on the operating rotor. However, 
there is no method for measuring'the latter derivative* directly. 'Rather,- each 
of the independent variables must be perturbed and the corresponding changes 
in the controlled variables noted. Furthermore, the independent variable must 
be perturbed in such a way that the changes in the controlled variable caused 
by these perturbations can be distinguished from other normally occuring 
changes in the controlled variable. This can be accomplished by using a per- 
turbation of sufficient amplitude or by accomplishing the perturbation with a 
recognizable modulation. In either case, a sufficient time must be provided 
for the dynamics of the rotor system to respond to the perturbation. 

The sensitivity of the optimization function to changes in the independent 
variable can be found in one of two ways: 

ML . 1 r ^ (v - V ) 

Ax, x ,2 - k2 kl 


where: 

3P 

is the derivative of the optimization function with respect 
'’"k to yj,, evaluated at the current operating point p) 
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or: 


AP 

Ax. 


^•2 


“ ^*1 


f(yi2> y22’ ^32* ^42^ - f(y-|p )/2V ^31’ 
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>] 


where: 

f(yi,*s) is the optimization function evaluated at the current opera- 
^ ting point (1) and the perturbation point (2), 

Theoretically, both methods are the same; however, there are some practical 
differences. For instance, the accuracy and resolution requirements on the 
calculation of the second are more stringent than on the first because the 
second involves the difference of two large numbers. However, depending on 
the nature of the optimization function, the storage of coefficients and the 
method of calculation may be simpler for the second than for the first. With 
either method, the resolution and accuracy requirements on the measurements of 
the controlled variables are the same. 

One method of perturbing the independent variables is to modulate them with a 
frequency, low with respect to the rotor dynamics, and then observe the con- 
trolled variables for a sufficient number of modulation cycles to separate the 
changes caused by the perturbation from the other changes. In general, the 
number of cycles required is inversely proportional to the amplitude of the 
perturbation. (The limiting case is 1/2 cycle; i.e., the independent variable 
is deflected from its value at the operating point for a sufficient time for 
the rotor dynamics to respond.) Since the rotor dynamics are fixed, the time 
to make a measurement is, in effect, inversely proportional to the amplitude 
of the perturbation. If some form of modulation is used, it would be possible 
to perturb several of the independent variables simultaneously and then separ- 
ate the responses of the controlled variables into the proper components cor- 
responding to each of the perturbed independent variables. However, if several 
are perturbed simultaneously, the amplitudes will combine and, therefore, the 
maximum perturbation amplitude allowed on each independent variable is inver- 
sely proportional to the number of variables that are perturbed simultaneously. 
Therefore, simultaneous perturbations will not decrease the time required to 
accomplish the measurements. 

The discussion above indicates that several rotor time constants are required 
to acquire the full set of derivatives. In addition, time must be allowed for 
the effect of the increment of the selected independent variable to take place. 
Therefore, because of its requirement to have all six derivatives, the maximum 
slope search described above may not be' the minimum time search. 

3.5.3 Sequential Step Search . The maximum slope search spends six time inter- 
vals determining derivatives before incrementing an independent variable 
towards the optimum point. During some of those six intervals, the system 
would have been incremented in such a way as to move closer to the "best" 
point. Therefore, the search time can be reduced if this possibility is recog- 
nized. This can be accomplished by noting whether the optimization function 
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decreases or increases as a result of each increment of an independent vari- 
able. The increment is retained if the optimization function decreases and 
remdved if it increases. The derivatives can also be determined when these 
increments are made. There are a variety of search techniques using this con- 
cept. They differ in the basic sequence used to increment the independent 
variables and in the logic used to modify that sequence in response to the 
results obtained from the increments. Possible criteria for establishing a 
basic sequence are: 

a. Increment each independent variable in turn, i.e., all with the 
same frequency 

b. Increment the most important independent variables more frequently 

. c. Increment the independent variable with the largest derivative 
most, frequently 

Inspection of the coefficients in Table 13 indicates that there is no apparent 
basis for establishing certain independent variables as being more important 
than others. Therefore, b. has not been considered. An approach based on c. 
cannot be accomplished without knowledge of the derivatives. These are con- 
stantly changing as the operating point moves and, therefore, must be deter- 
mined at some minimum interval. This can be accomplished with a method based 
on a. The method chosen for the purposes of this feasibility study uses a com- 
bination of a and c. A basic 12-interval cycle is established. During the odd 
numbered intervals, each of the six independent variables is incremented in 
turn. The direction of the increment is determined on- the basis of the last 
slope determined for that independent variable. The even numbered intervals 
are used to increment the independent variable with the largest derivative. 

A derivative is determined in each interval and stored for future decisions. 

Variations of this strategy are also categorized by the action taken when it 
is determined whether or not the increment of .an independent variable 
increases or decreases the optimization function. Some of the possibilities ' 
are: • ' . • 

. I * ’ 

a. Remove the increment if the optimization function increases 

b. Retain 1/2 the increment if the optimization function increases 
and if the slope changes sign, the assumption being that the 
increment caused the operation to move past the minimum point 

c. Try an increment of the opposite sign if the optimization func- 
tion increases and the slope has changed sign. 

\ 

The effects of these variations and others have been studied in the simulation 
studies described in Section 4. 

3.5.4 Magnitude of the Increment . There are two considerations that call for 
a large increment of the independent variables. These are better immunity to 
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noise and other random variations in the controlled parameters, and more rapid 
response of the system to changes in operating condition of the helicopter. 

The considerations for a smaller increment are a closer approach to the minimum 
point and less likelihood that the effect of the perturbations become objec- 
tionable to the operation of the helicopter. Some of each of these consider- 
ations can be accommodated by making the increment variable. The advantage of 
this complexity is dependent upon the noise levels that will be encountered, 
rate of change required of the multi cyclic control components to accommodate 
changes in helicopter operation, and the sensitivity of the search to the size 
of the increment. Some of these effects are Investigated in the simulation 
studies described in paragraph 4.2. However, the determination of the magni- 
tude of the increment will be largely dependent upon experience with an oper- 
ating rotor. 

3.5.5 Choice of Independent Variables . While it was stated in paragraph 3.1 
that the independent variables would be the sine and cosine components of the 
2, 3 and 4/revolution pitch angle, consideration should also be given to using 
the magnitude and phase angle of the 2, 3 and 4/revolution pitch angle as the 
independent variables. The latter approach might reduce the search time in 
response to changing operating conditions. This would be true if the best 
phase angles of the multicyclic control angles are a function of the direction 
of the rotor thrust vector and if the best magnitude of the multicyclic con- 
trol angles are a function of the magnitude of the rotor thrust vector. In 
order to determine the advantages of this approach, rotor characteristics 
would be required describing a variety of operating conditions. 

There are also disadvantages to using magnitude and phase angle as the inde- 
pendent variables. For instance, when the magnitude of one harmonic is near 
zero, changes in the phase of that harmonic will have little or no effect. 
Furthermore, at a given phase angle, the best magnitude may be near zero. 

These two conditions can cause a suboptimum "depression" in the optimization 
function surface. The magnitude of this problem is dependent upon the rotor 
characteristics for various operating conditions and the characteristics of 
the optimization function. 

In view of the potential problems and the lack of any clear indication of a 
specific advantage to using magnitude and phase angle as independent variables, 
the sine and cosine components are used as the independent variables in this 
feasibility study. 

4.0 SIMULATION STUDIES 

Many of the factors that must be investigated in the course of developing a 
preliminary design for the feedback control system cannot readily be treated 
analytically. Therefore, two computer models have been developed to allow 
these effects to be investigated. The first of these is based on the maximum 
slope method for selecting the operating point. Since this represents the 
most effective search procedure {without regard to the determination of the 
derivatives of the rotor characteristics), it is used to investigate those 
effects that are not associated with a particular search method. It is also 
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used as the standard against which to compare the results of other search 
methods. The model is described in Appendix C, The other model allows the 
investigation of particular search strategies and the way in which they are 
affected by various system parameters. This model is discussed in Appendix D. 
The results., using each of these models, are discussed below. Unless noted 
otherwise, each model used the optimization function of Figure 3 and an incre- 
ment of 0.5°. 


4.1 Effect of Various Optimization Functions 

The relative value of the various controlled parameters, as reflected by the 
optimization function, determines the operating point selected by the control 
system action. The effect of the optimization function on this selection pro- 
cess is dependent upon the rotor characteristics. The interaction of various 
optimization functions with the given set of rotor characteristics has been 
investigated using the maximum slope model. This model was chosen so that the 
results will be the least dependent on the particular search strategy used. 

• 

Two basic optimization functions were used in this investigation. The first 
places equal value on each of the controlled parameters. This is represented 
by the following equation: 


P = 100 


f BM ^ HP . “max . RV ] 
7000 750 12 300j 


The second is the optimization function presented in paragraph 3.2 and des- 
cribed in Figure 3. The variations investigated include various combinations 
of the terms from these two functions, as well as changes in the constants 
used in the second. The results are judged in terms of the operating point 
selected by each function as described by the values of the four controlled 
variables. The results are summarized on Table 14. The operating point with 
no multicyclic control has been included as a basis for comparison with the 
operating point selected by the various optimization functions. 


The first case is the linear optimization function where each of the con- 
trolled variables, normalized to a particular value, has equal weight. This 
shows that a significant reduction in the hub shear (RV) can be achieved with 
only a small increase in the other three controlled variables. 

Cases 2, 3 and 4 show the cummulative effect of successive changes of the 
terms associated with hub shear, maximum angle of attack and horsepower, from 
those of Case 1 to the optimization function described in paragraph 3.1. The 
change in the hub shear term places a greater importance on hub shear with the 
result that it is decreased at the expense of increases in the other three 
parameters. The change in the angle of attack term eliminates the penalty of 
angles near 12° so that the angle of attack is further increased, providing an 
additional reduction in hub shear. The change of the horsepower term 
increases the penalty of 771 HP. Therefore, horsepower is decreased and hub 
shear increases. The change in maximum angle of attack is just the conse- 
quence of the different selected operating point. The change in the bending 
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TABLE 14, EFFECT OF VARIOUS OPTIMIZATION FUNCTIONS 


SELECTED OPERATING POINT 


INDEPENDENT VARIABLES CONTROLLED VARIABLES 


OPTIMIZATION FUNCTION 

2S 


2C 

3S 

3C 


4S 

4C 

BM 

HP 

“m 

RV 

Base Line: No Multicyclic Pitch 

0 


0 

0 

0 


0 

0 

4411 

745 

IKl 

297 

1, 

D - BM .HP . “m , RV 
^ ■ 7000 750 12 500 

2 

- 

.5 

- .5 

0 

- 

1 

.5 

4549 

759 

1K7 

211 

2. 

Change RV term to that of Figure 3 

2.5 

- 

.5 

.5 

0 

- 

K5 

1 

4653 

765 

12.0 

197 

3. 

Also change term to that of Figure 3 

3 

- 

.5 

.5 

0 , 

- 

1.5 

1 

4723 

771 

12.4 

186 

4 

Also change HP term to that of Figure 3 

2.5 


.5 

0 

0 

- 

K5 

1 

4653 

765 

12.0 

197 

5. 

Optimization Function of Figure 3 
(i.e., also change BM to that of Figure 3) 

2.5 

- 

.6 

.5 

0 

- 

1.5 

1 

4653 

765 

12.0 

197 

6. 

#5 without linear BM term 

2.5 

- 

.5 

0 

0 

- 

K5 

1 

4683 

766 

12.1 

195 

7. 

ilfS with HP threshold increased from 750 to 800 

3 

- 

.5 

0 

0 

- 

1.5 

1 

4723 

771 

12.4 

186 

8. 

#5 with threshold increased from 12 to 13 

3 

- 

.5 

.5 

0 

- 

1.5 

1 

4688 

770 

12.3 

188 

9. 

Kb with HP^ = 800, ^ 

4.5 

- 

1 

0 

0 

- 

1.5 

1 

5054 

792 

13.6 

155 

TO. 

KS without linear BM term, HP^ = 800 

3.5 

- 

1 

0 

0 


K5 

1 

4880 

777 

12.7 

173 

IK 

Kb without linear BM term, = 13 


- 

1 

0 

.5 

- 

1.5 

1 

4794 

770 

12.4 

180 ‘ 

12, 

Kb without linear BM term, HP. = 800, a ^ = 13 

4.5 

- 

1.5 

0 

0 


K5 

1 

5165 

792 

13.5 

153 
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moment term results, in Case 5, in the full optimization function of Figure 3. 
Since this new term has essentially the same value as the original term at a 
bending moment of 4653 in-lb, the operating point is unchanged. 

Cases 6, 7 and 8 show the effects of making changes in the bending moment, 
horsepower and maximum angle of attack terms. In Case 6, the linear portion 
of the bending moment term is removed. The operating point then selected indi- 
cates that, with these rotor characteristics, allowing greater bending moments 
does not significantly decrease the hub shear that can be obtained. However, 
Cases 7 and 8 indicate that increasing the threshold for either horsepower or 
maximum angle of attack does allow a decrease in hub shear to be obtained. 

The fact that the horsepower and maximum angle of attack are essentially the 
same in Cases 7 and 8 is .caused by the particular rotor characteristics being 
used and not the optimization function. The optimization function is assigning 
an essentially equal penalty in either case whether the penalty is caused by 
horsepower being above its threshold or angle of attack being above its thresh- 
old. 

Cases 9, 10 and 11 show these three changes being applied in pairs. In Case 9, 
the thresholds for both horsepower and angle of attack have been increased and, 
as would be expected, another significant decrease in hub shear is achieved. 
Cases 10 and 11 again show that bending moment does not represent a serious 
constraint with these rotor characteristics. Although the elimination of that 
term does allow a greater decrease in hub shear than the increase of the two 
thresholds by themselves. 

Case 12 incorporates all three of the changes described above and, in effect, 
represents the addition of the bending moment deletion to Case 9. As indi- 
cated, this does not provide much further reduction in hub shear. 

In conclusion, it can be seen that with these rotor characteristics, signifi- 
cant decreases in hub shear, and the resulting body accelerations, can only be 

achieved at the price of increases in horsepower, maximum angle of attack and, 
to some extent, bending moment. Therefore, the optimization function selected 
must accurately reflect the relative value of these parameters for the multi - 
cyclic control system to be effective. 

4.2 Magnitude of the Increment 

The factors to be considered in chosing the amplitude with which to increment 
the various independent variables were discussed in paragraph 3.5.4. The 
effects of various increment magnitudes have been investigated with the maxi- 
mum slope model in the absence of other factors, such as noise and the peculi- 
arities of particular search strategies. Some of the results are shown on 
Figure 4. The optimization function of Figure 3 was used and the model was 
started at the operating point corresponding to no multicyclic control. The 
increment was first set at 1° and the model allowed to run until it reached 
its best operating point. The increment was then changed to .5® and the model 
was allowed to continue until it again reached its best operating point. Two 
additional cycles were taken with the increment at .2'* and .1°, respectively. 
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Figure 4. Effect of Increment Size. 
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These results indicate that, with the optimization function and rotor charac-- 
teristic used and without considering noise in measurement of the rotor char- 
acteristics, there is little benefit in decreasing the increment below .5°. 

4.3 Effect of Search Strategies 

The general considerations for search strategies were discussed in paragraph 
3.5, with particular emphasis on sequential step searches in paragraph 3.5.3. 

The simulation model described in Appendix D was set up to investigate the 
characteristics of the various alternatives for sequential step searches. Two 
sequences were investigated: one with a basic cycle of six intervals,’ wherein 

each of the independent variables was incremented, in turn; and the other with 
a basic cycle of twelve intervals, wherein the six independent variables were 
incremented in the odd numbered intervals and the independent variable having 
the maximum slope was incremented in the even numbered intervals. Within each 
of these sequences, four alternatives on the action to be taken if the optimi- 
zation parameter increases as a result of an increment were investigated. 

These four alternatives are: 

A. Retain the increment and move to the next interval in the cycle 

B. Remove the increment and redetermine the operating point without 
the increment 

C. Remove one-half the increment and take the new value of the opti- 
mization parameter to be the average of the values with and 
without the full increment 

D. Try an increment of the opposite sign. If the optimization 
parameter still increases, remove the increment and redetermine 
the operating point without the increment. 

In each case, the slope (AP/AX) is determined in each interval and stored for 
use in determining the direction of the Increment the next time the corres- 
ponding independent variable is to be incremented. In addition, in method D, 
if the increments with opposite signs both cause the optimization parameter to 
increase, the slope is set equal to zero, with the sign equal to that of the 
smaller of the two slopes. 

In methods B and D, the redetermination of the operating point with the incre- 
ment removed is necessary in order to allow the system to follow changes in 
rotor trim conditions. It is also necessary to account for the effects of 
noise and other random variations. 

The results are summarized on Table 15. This shows, for each of the cases, 
the number of iterations for the optimization parameter, P, to reach 60, 50 
and its minimum. The minimum value of P’ is also given. The remaining four 
values in the table are the average values and the standard deviation of the 
optimization parameter, P, and the hub shear, R^, after steady state has been 

reached. The particular value of R^ is, of course, a function of the 
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lA Continue 
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46 
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19 

32 
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49 
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optimization parameter- However, the value of the standard deviation of R^ is 

an indication of the stability of the operating point achieved. A plot of the 
operation of method 2D is shown on Figure 5. 

Comparison of the four methods, wherein each independent variable is incre- 
mented in sequence, shows that method A has the largest average optimization 
parameter and the largest standard deviations.- This is caused by the incre- 
ments that drive the system away from the "best" operating point remaining 
in existence for the six intervals of the basic cycle. In case B, where the 
increment is removed if P increases, the steady state operation is improved 
without any significant change in the time to reach the steady state condition. 
Removing one-half the increment, as in method C, provides a slight reduction 
in the time to reach steady state. However, the steady state condition is more 
like that of method A. Method D was an attempt to reduce the time to reach 
steady state while retaining the steady state characteristics of method B. 
However, the improvement in time was not achieved. 

The sequence wherein the sequential incrementing of the independent variables 
is alternated with the independent variable providing the maximum slope, gen- 
erally reduces the number of iterations necessary to reach the steady state 
condition. Method D showed the greatest reduction and does provide a signifi- 
cant decrease in the number of iterations to achieve steady state while 
retaining the steady state characteristics of method B. The improvement in 
the steady state conditions for methods 2A and 2C over lA and 1C, respectively, 
is probably due to removing the bad excursions from the "best" operating point 
sooner through the provisions for incrementing the independent variable pro- 
viding the maximum slope. 

The effect of increment size is shown by case 2D', which is the same as 2D, 
except that the increment is 1 degree (the increment is 1/2 degree for all 
other cases). While it might have been expected that the time to reach steady 
state would have been reduced, such was not the case. Undoubtedly the results 
would have been different if a greater amount of multicyclic control were 
required to reach the "best" operating point. Whether or not this would ever 
be an advantage is dependent on how fast the operating point must be changed 
in response to rotor trim changes. The deterioration of the steady state con- 
dition is caused by the larger excursions from the "best" operating point 
caused by the larger increments. 

Methods 2C and 2D appear best, with 2C providing somewhat faster initial move- 
ment towards the "best" operating point, but 2D providing slightly better 
steady state conditions. 

The final choice of a search strategy is dependent upon the requirements for 
the system to follow changes in rotor trim conditions. This is, in turn, 
dependent upon the rotor characteristics under a variety of trim conditions 
and the characteristics of the final optimization function. 
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Figure 5. Response to Search Strategy 2D. 
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4.4 Noise Investigation 

The control system must make measurements of the dependent variables to provide 
feedback to the control process. In addition, the measurements of these depen- 
dent variables also provide the basis for determining the derivatives of the 
optimization function with respect to each of the independent variables. How- 
ever, the instrumentation used to measure these dependent variables is sub- 
jected to various external effects. In addition, the dependent variables them- 
selves are perturbed by sources not related to the increments of the indepen- 
dent variables. Both of these effects constitute noise which can mislead the. 
control system. 

To investigate these effects, a noise- generator was included in the simulation 
model described in Appendix D. A random process is used to- generate noise with 
an amplitude that is normally distributed about zero. These noise amplitudes 
are then added to the calculated values of the dependent variables. The noise 
generator has been arranged so that the RMS amplitude of the noise added to a 
particular dependent variable is equal to a given percentage of the threshold 
value used in the optimization function for that variable. 

Runs of the simulation model were made with RMS noise amplitudes of .2%, . 5 %, . 
1%, 2% and 5 % applied to all of the dependent variables. The results are shown 
on Table 16. No results are shown for the 5% case, since the system did not 
show signs of reaching a stable operating condition during the first 175 iter- 
ations. P was calculated before and after the noise was added to the dependent 
variables, and information about both is provided on Table 16. The R^ infor- 
mation is also given without noise, since that represents the actual R^. 

Since the noise generation is a random process, many runs at each of the opera- 
ting conditions would be necessary in order to draw definitive conclusions. 

Such has not been done within the scope of this feasibility study. However, 
from Table 16, it would appear that noise begins to have an effect at 1% RMS. 

The 2 % case was interesting in that the system seemed to find one operating 
point after approximately 65 iterations. This point was held for 30 or 40 
iterations and then the system moved to a different operating point, as indi- 
cated on Table 16. The run was terminated after 200 iterations, and at that 
point there was an indication that the operating point was again shifting. 

A detailed investigation of some of the data produced by the simulation model 
indicates that the noise added to the horsepower is a major contributor to the 
instability evidenced by the model. A review of the derivatives of all of the 
dependent variables with respect to the -independent variables also indicates 
that' horsepower is most likely to be a problem. These derivatives are given 
on Table 17 for the initial point, that is the operatihg point with no multi- 
cyclic control, and the "best" point. These derivatives are calculated from 
the data given in Table 13. This data shows that, except for the sine component 
of the 2/rev at the "best" point, the derivatives of horsepower are very small. 
Changes of the multicyclic components of one-half degree provide only a very 
few tenths percent change in horsepower. Thus, with noise of ^% RMS and 
"signal" of only a few tenths of a percent, the signal to noise ratio is very 
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TABLE 16. EFFECTS OF NOISE 


RMS NOISE 

% 

ITERATIONS TO 
STEADY STATE 

BEFORE NOISE 
P ap 

WITH 

P 

NOISE 

^P 

BEFORE 

''v 

NOISE 

■‘^Rv 

0 

40 

48 

1.8 

48 

1.8 

189 

6 

.2 

50 

49 

2.0 

49 

2.4 

188 

6 

.5 

50 

49 

1.5 

49 

1.7 

198 

5 

1.0 

55 

51 

3.3 

55 

9.6 

204 

7 

2.0 

65 

66 

4.3 

69 

7.4 

241 

9 


135 

52 

3.5 

55 

9.2 

204 

9 

5.0 

Unstable 
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TABLE 17. VALUES OF DERIVATIVES 



1 \ 

INDEPENDENT VARIABLES 
COMPONENT VALUE 

DERIVATIVES OF DEPENDENT VARIABLES 
WITH RESPECT TO INDEPENDENT VARIABLES 

%ax 

(% OF THRESHOLD/DEGREE) 

Initial Point 






2s ' 

0° 

- 2.0 

0.2 

1.6 

5.9 

2c 

0° 

- 2.9 

- 0.2 

5.0 

3.1 

3s 

0° 

0.5 

- 0.4 

- 0.6 

9.0 

3c 

0“ 

0.8 

- 0.2 

- 0.1 

0.5 

4s 

0° 

- 0.1 

0.0 

- 0.3 

16.3 • 

4c 

0° 

0.8 

0.1 

- 1.1 

5.3 

“Best" Point 






2s 

3.0° 

1.4 

1.4 

6.5 

2.5 

2c 

1 

0 

O 

o 

- 3.6 

- 0.1 

1.7 

0.6 

3s 

- 0.5° 

- 1.5 

- 0.4 

- 1.8 

1.0 

3c 

0.5° 

0.1 

0.0 

2.3 

2.4 

4s 

- 1.0° 

- 2.4 

- 0.6 

- 2.5 

2.0 

4c 

0.5° 

- 0.7 

0.4 

- 0.1 

1.8 

1 ^ -- - ^ 
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low. This is aggravated by the fact that at the "best" operating point the 
derivative of the optimization function with respect to horsepower is very 
Targe. Therefore, the system sees large changes, in the optimization parameter 
and; attributes these to the incremental changes being made to the independent 
variables when, in fact, the Targe changes in optimization parameter are a 
resu.1t of the noise on the horsepower signal. That the system functions at 
all in the face of 1% noise can be attributed to the fact that the system con- 
tinuously makes many, many measurements of the dependent variables and so the 
effects of noise tend to average out, at least to some extent. 

The results above indicate that .noise .must be given special consideration 
during the detail design. • The effects of noise will be influenced by the 
rotor characteristics and the characteristics of the optimization function 
finally chosen; These characteristics and the characteristics of the noise 
anticipated in the measurement of each of the dependent variables will deter- 
mine what type of signal processing will be necessary to reduce noise to 
acceptable levels. 

4.5 Dynamic Considerations 

In feedback control systems, action is taken on the basis of response to the 
previous action. Therefore, the total system characteristics are dependent 
upon the dynamics of the response of the controlled element to variations in 
the independent variables. Generally, the dynamic response of the controlled 
element is given and cannot be considered a design variable to the designer 
of the control system. In this multicyclic control system, the major design 
variable that influences the total system dynamics is the rate at which the 
independent variables are incremented, or the iteration interval. 

In order to investigate the effects of the iteration interval, a way is needed 
to characterize the dynamic response of the rotor system and the instrumenta- 
tion and signal processing used in the measurement of the dependent variables. 
For the purposes of this feasibility study, this response has been character- 
ized as a first order time lag. This is adequate for establishing the gross 
effects of various iteration intervals. It certainly would not be adequate 
for an investigation as to how the dynamic characteristics of the control sys- 
tem might be adjusted to compensate for the dynamics of the rotor. 

The effects of dynamic response have been incorporated in the simulation model 
by adjusting the calculated values of the dependent variables in accordance 
with the conventional equation for the response of the first order system as 
shown below: 

y(t) = y,, - (y^ - y„) e'*/'' 
where: 


47 



iimxM COR^PORATION 

OLD WINDSOR ROAD, BLOOMFIELD. CONNECTICUT C6003 


A 


CORPORATION 


Report No. R-1494 
21 January 1977 


= calculated value of y 
y^ = value of y at beginning of interval 
.T = time constant equivalent of the dynamic response 
t = length of iteration interval 

This adjustment provides the value of the dependent variables at the end of 
the' iteration cycle and it is these values that are used to calculate the opti- 
mization function which, in turn, is used to determine the slopes and the sub- 
sequent control action. While the model has the capability of incorporating a 
different time constant for.each of the dependent variables, the work in this 
feasibility' study has been' confined to the use of a single time constant. It 
will be noted that the response of the system can be characterized by the 
ratio of the iteration interval to the rotor time constant. 

Runs were made with the simulation .model of Appendix D under conditions ranging 
from an iteration interval very long with respect to the rotor time constant to 
one which was 10% of the rotor time constant. The results are summarized on 
Table 18. The important parameters are the time to reach a more or less steady 
state and the conditions that exist during that steady state. The time para- 
meter on Table 18 is expressed in terms of the rotor time constant and is 
obtained from the product of the number of iterations and the ratio of the 
iteration interval to the rotor time constant. Conditions during steady state 
are characterized by the mean and standard deviation of the -optimization para- 
meter, P, and the hub shear, R^. Time, the mean value of P, and the ratio of 

the standard deviation to the mean of hub shear are also plotted on Figure 6 
as a function of the ratio of iteration interval to rotor time constant. 

It will be noted that the time to reach steady state continues to .decrease as 
the iteration interval is decreased, until the interval become's less than 15% 
of the rotor time constant. Below that point, the system did not arrive at a 
steady state condition within 300 iterations, or 30 rotor time constants’. The 
scatter of the points on Figure 6 is exactly reproducible for identical condi- 
tions and is, undoubtedly, caused by the particular sequence of events that 
take place as a result of a particular set of conditions. More consistent 
results would probably be obtained on the basis of averages derived from using 
a variety of operating points. An example of the action that can take place 
in response to a particular set of circumstances occurred when the ratio of the 
iteration interval to the rotor time constant was .24. In this case, the oper- 
ation settled down to the conditions shown on Table 18 after approximately 
126 iterations. After an additional 325 iterations, a new operating point was 
found, wherein the mean value of the optimization parameter was 52, and its 
standard deviation 1. 

Figure 6 indicates that the performance of the system begins to degrade as the 
iteration interval becomes about one-half the rotor time constant. However, 
the final determination of this threshold must await further definition of the 
rotor operation under other rotor trim conditions, and determination of the 
dynamic response associated with each of the measured variables. 
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1 - 


• ITERATIONS 

t/x 

P . 

min 

TO P . 
rtnn 


1 

46 

52 

3 

.95 

46 

43 

2 

.86 

46 

44 

1.5 

.78 

47 

45 

1.2 

.70 

46 

54 

1.0 

.63 

48 

48 

.8 

.55 

49 

40 

.6 

.45 

47 

41 

.4 

.34 

48 

39 

.3 

.26 

■ 52 

56 

,24 

.21 

56 

42 

.2 

.18 

59 

96 

.15 

.14 

60 

44 

.1 

.095 

60 

118 

*In 

units of rotoi 

' time 

constants. 
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4.6 Combined Noise and Dynamic Effects 

In paragraphs 4.4 and 4.5, it was indicated that either noise or dynamic 
effects can influence the stability of the feedback control system and the 
quality of the operating point. These effects became noticeable with an RMS 
noise of 1%., or an iteration interval of about 1/2 of the effective system 
time constant. Therefore, these two effects have been tried in combination to 
determine the extent to which they interact. 

The simulation model of Appendix D was run with an RMS noise of 1 % and an iter- 
ation interval of 1/2 the equivalent rotor time constant. While the system 
did move towards the previous operating point for 60 to 80 iterations, it then 
began to move rapidly away from that operating point. Reducing the noise by a 
factor of 2 had very Tittle effect. Doubling the iteration interval did pro- 
vide an improvement in that the system did come close to the previous operating 
point and remained close to it for several hundred iterations. However, it 
then began moving away rapidly. 

Several other combinations of noise and iteration- interval were tested, and the 
results are summarized on Table 19. This shows that there is considerable rein- 
forcement of the effects of noise and iteration interval with, perhaps, the 
system being somewhat more sensitive to the iteration interval. It also indi- 
cates that there must be a significant increase in the iteration interval or 
decrease in the noise level over their respective "threshold values", as deter- 
mined independently, for acceptable performance. A combined factor of at least 
4, and more likely 8 or 16, appears necessary. 

5.0 PRELIMINARY DESIGN • 

5.1 Impl emen tati on Concepts 

The basic algorithm for multi cyclic control and its performance in conjunction 
with the rotor characteristics of Table 11 and the optimization function of 
Figure 3 has been described in the previous sections. This section is con- 
cerned with the implementation of this control algorithm and its interface 
with the existing MFS/HPR Module System. 

The calculations and decision points involved in the control algorithm strongly 
suggest the use of digital concepts in its implementation. Furthermore, since 
flexibility will be an important attribute for the test configuration, program- 
mable digital concepts should be used. This can be accomplished with either a 
mini- or a microprocessor, depending upon availability and the degree of flex- 
ibility determined necessary. A teletype or CRT terminal in conjunction with 
the .processor can provide the necessary operator inter'face, although other 
methods can also be considered. This provides flexibility in the displays and 
controls available to the operator through changes in the processor program. 

There is an analog control system in existence with provisions for generation 
of the second, third and fourth harmonics of the rotor azimuth angle, provision 
for modulating, or changing the amplitude and phase angle, of these harmonics 
with manual inputs, and provisions for summing the' modulated harmonics to 
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TABLE 19. COMBINED NOISE AND DYNAMIC EFFECT 

(ITERATION INTERVAL) 
(EFFECTIVE TIME CONSTANT) 

RMS NOISE .5 1 2 


1 % 

Unstable 

Divergent 

Fair 

1/2 % 

Unstable 

Marginal 

Fair 

1/4 % 

Marginal 

Fair 

Good 


Categories : 

Unstable 

Divergent 

Marginal 

Fair 

Good 


Does not reach or does not hold a stable 
operating point. 

Reaches and holds a stable operating point for 
a few hundred iterations and then diverges. 

Reaches and holds a stable operating point, but 
wanders away from it for short periods 

Reaches and holds a stable operating point, 
such that the optimization parameter stays 
within 20% of the "value" without noise and 
dynami cs . 

Operation essentially uneffected by noise and 
dynamics. 
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generate the specific control signals for each of the four blades. These con- 
trol signals, in turn, drive the blade servos associated with each blade. This 
system has shown itself to be more than adequate for the puposes of the test 
configuration. Therefore, the feedback control system should interface the 
existing system at the modulation step. The control processor can provide mod- 
ulation signals for each of the harmonics that replace the function of the 
existing manual inputs. However, it is more convenient to provide the modula- 
tion signals and implement the modulation process in terms of the sine and 
cosine components of each harmonic, rather than the amplitude and phase, as is 
done with manual control. 

A block diagram of the preliminary design concept is shown in Figure 7. Sig- 
nals from each of the. sensors are conditioned by analog circuitry for input 
to the control processor through an analog to digital converter. {When the 
total signal processing task has been established on the basis of more complete 
knowledge of the characteristics of each of the measured signals, a determina- 
tion will be made as to which parts can better be done in the analog signal 
conditioning circuits and which parts might better be done digitally within 
the control processor.) The signal conditioning for the body acceleration 
signals is the most complex, in that the requirement is to extract the ampli- 
tude of the fourth harmonic of a weighted average of these signals. The fourth 
harmonic can be extracted by a synchronous detector using the sine and cosine 
components of the fourth harmonic as generated by the existing harmonic gener- 
ator. It may also be desirable to extract a particular harmonic from the con- 
trol load signal which is representative of bending moment. The average of 
the horsepower signal and the peak of the pitch link load signal is determined. 

The processor senses each of the measured inputs and implements the control 
algorithm to generate the amplitude of the sine and cosine components of each 
of the three harmonics. These six signals are then used, to modulate the har- 
monics of rotor azimuth as developed by -the existing harmonic generator and, 
thus provide signals that replace those generated by the present, manually 
controlled modulator. These signals are then summed in the existing ampli- 
fiers to provide the specific control signal for each of the four blades. 

The processor also provides the interface with the operator through a teletype 
or CRT terminal. Other methods may also be considered when the complete 
requirement for operator interface is established. Among requirements to be 
considered are display of the measured variables and the calculated indepen- 
dent variables, initiation and termination of feedback control, change of con- 
trol parameters such as iteration interval and increment size, and override of 
feedback control with manual control. 

Other degrees of flexibility which should be provided^ but not through the 
operator interface, include changes in the signal processing characteristics 
and the control algorithm. 

5.2 System Test Sensors 

As stated in paragraph 3.1, the parameters selected for control in the feed- 
back system are Bending Moment, Horsepower, Maximum Angle of Attack and Hub 
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Shear. All studies to date have indicated that these parameters should be the 
primary controlling factors indicating satisfactory multi cycli.c performance. 

It is possible" that data from wind tunnel tests may indicate that additional 
parameters may be sensed to limit the harmonic inputs to maintain all systems 
within acceptable load levels. The feedback system, which is the subject of 
this study, does not consider these additional limiting parameters; however, 
they could be added when necessary. 

Direct measurement of the parameters listed above presents some problems. 
Ultimately, a flight version of the feedback system must be considered which 
makes the measurement of bending moment Impractical from a reliability stand- 
point. Horsepower, angle of attack and hub shear cannot be measured directly 
in either an aircraft or the wind tunnel model. It is necessary, for all of 
the above parameters, to select control parameters that are either directly or 
indirectly related to the desj red parameters. . 

Blade bending moment can only be measured directly by strain gage bridges. 

While this would be acceptable for a wind tunnel model, the installation of 
strain gage bridges on flight blades that may be exposed to varying weather 
conditions is not satisfactory for a reliable control system. There is a high 
probability that blade bending moments are related to control loads. The rela- 
tionship must be verified and defined through analysis of wind tunnel data. 

If there is, in fact, a relationship, reliable load cells can be installed in 
the fixed control system. For wind tunnel test, it is more practical for cast 
considerations to use existing instrumentation that measures pitch link load. 

Horsepower is a more straightforward parameter with which to work. In the 
wind tunnel model, rotor shaft torque can be used, whereas in a helicopter 
transmission, torque pressure, for which standard sensors are available, would 
be used. 

Blade angle of attack cannot be directly measured, but there is a relationship 
between angle of attack and pitch link load. Again, in a flight helicopter, 
it would be more desirable to work with load cells in the non-rotating control 
system. 

Hub shear has always presented a measurement problem. In the case of multi - 
cyclic control, the interest in vibratory hub shear reduction is due to the 
resultant reduction in fuselage vibration. Vibration measurement by accelero- 
meters will provide the initial control measurement. It is uncertain, due to 
data availability, whether or not the accelerometer signals can be fed directly 
to the feedback circuitry or if the vibration levels must be first related to 
hub shear through methods of Force Determination. Tbe number and location of 
the accelerometers is also uncertain until more data becomes available. For 
the purpose of this study, the accelerometer signals are fed directly to the 
feedback system. 

5.3 Processor Functions 

A flow diagram of the processor functions associated with the control algorithm 
is shown on Figure 8. (The other major function of the processor, which is not 
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SET A, B, C = 1 
INITIALIZE X's 
INITIALIZE S's 
SET MAGNITUDE I 


CALCULATE 
OPT. PAR., PI 


A = 1 


MEASURE 

INPUTS 



1 





SET B = 1 

PROCESS 

INPUTS 

! 

FIND N FOR MAX S 


*SET SIGN OF TO SIGN OF 


THE SMALLER OF S. AND (P2 - Pl)/I 



S . <.01 


SET A 2, P2 = PI 
SET SIGN OF I = SIGN OF - S, 

SET + I 


Figure 8. Processor Control Functi 
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shown, is the servicing of the operator interface.) Three branches are shown. 
The "Initiate" branch is entered each time the feedback control function is 
initiated. Its purpose is to provide the initial values of the six components 
of multicyclic control (X's) and the initial values of the slopes of the opti- 
mization function with respect .to each of these components (S‘s). The latter 
wilT usually be initialized at zero. Two flags, A and B, and a counter, C, 
that are used to control the operation of the control algorithm are initialized 
to the 1 state. The magnitude of the increment to be' used in conjunction with 
the independent variables must also be set. 

The "Input" branch is associated with measuring each of the inputs and accom- 
plishing that portion of the signal processing that has been assigned to the 
control processor. Depending upon the nature of the signal processing to be 
accomplished within the control processor, this branch may be entered more fre- 
quently than the control branch. 

The "Control" branch is concerned with the control function itself and is 
entered once each iteration interval. The first step is to calculate the opti- 
mization parameter from the measured and' processed inputs. Next, flag A is 
checked to determine if there is a valid basis for determining whether the 
optimization parameter has increased or decreased as a result of the most 
recently incremented independent variable. If there is, and if the optimiza- 
tion parameter has decreased, the current and previous values of the optimiza- 
tion parameter and the increment, I, are used to calculate the slope to be 
associated with the most recently incremented independent variable, Sj^. Flag 

B is then examined to determine whether the next independent variable to be 
incremented is to be selected on the basis of maximum slope or the numerical 
sequence. The counter, C, is used to keep track of the numerical sequence. 

The maximum slope criterion is not used if the maximum slope is less than .01. 
After setting the value of N (the independent variable to be incremented), the 
value of A and the baseline value of the optimization parameter, ? 2 i are set 

to the values appropriate for the next iteration cycle. The independent vari- 
able corresponding to N is then incremented in a direction dependent upon the 
sine of the slope associated with that independent variable. 

If, back in the beginning of the "Control" branch, the optimization parameter 
were found to have increased, and if flag A is not in the 3 state, steps are 
taken to increment the most recently incremented independent variable in the 
opposite direction. If this had already been accomplished (i.e., flag A in 
the 3 state), the last increment is removed and flag A is set to 1 so that a 
new baseline value of the optimization parameter is determined on the next 
iteration cycle. In addition, the sign of the slope associated with the most 
recently incremelited independent variable is set equal to the sign of the 
smaller of the slopes determined from the positive and negative increments. 

The magnitude of the slope is set near zero so that it is not considered in 
the determination of the independent variable with the maximum slope. 
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5.4 Modulator Function 

A schematic of one implementation of the modulator function is shown in Figure 
9. One set of inputs is the sine and cosine components of the second, third 
and fourth harmonics of the rotor azimuth angle generated by the existing har- 
monic generator. The second set of inputs is the required amplitudes of these 
components as determined by the control processor. The eight blocks marked M 
generate the product of their respective inputs. Pairs of these products are 
then summed to provide the signals required by the existing summing amplifiers 
for each of the four blades. The M blocks may be implemented in either of two 
ways. The first uses digital to analog converters on each of the six control 
processor outputs and eight, four-quadrant, analog multipliers. The second, 
uses eight ladder networks, each being switched by the appropriate digital 
output of the control processor. 

6.0 CONCLUDING REMARKS 

During this feasibility study, several of the concepts associated with multi- 
cycTic adaptive control have been investigated and a preliminary design for 
their implementation has also been generated. In the course of the investiga- 
tion, some questions have been uncovered which require further resolution; 
however, no fundamental problems were discovered. Th,e preliminary design also 
indicates that the concepts can be implemented with available technology. 

There are certain areas where further work is required to support the develop- 
ment of a detailed design. These areas may be grouped in four categories: 
the characteristics of the optimization function; the characteristics of the 
measured variables; control system stability; and rotor characteristics. 

6.1 Optimization Function Characteristics 

The simulation studies have shown that the optimization function not only 
affects the selected operating point, but also affects certain system operating 
characteristics. It is important to know the nature of the criteria (i.e., 
the tradeoffs among the controlled, or measured, variables) for describing 
the “best" rotor operation and the nature and numbers of controlled variables. 
Consideration should also be given to the necessity for changing these criteria 
as a function of the rotor trim conditions. 

.6.2 Characteristics of Measured Variables 

The characteristics of the optimization function will identify the kinds of 
variables that must be measured to provide the necessary feedback. Knowledge 
is required of the nature of the dynamic response of^each of these measured 
variables to changes in the independent variables, as well as the amplitude 
and frequency characteristics of the noise to be associated with each of 
these measured variables. The latter is particularly important in establishing 
the nature of the signal conditioning necessary to reduce the noise to accept- 
able levels. Since these characteristics are not peculiar to multicyclic con- 
trol, they may be developed on the basis of experience with other helicopter 
rotors and their instrumentation. 
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6. 3 Control System Stability 

Further study is required of the interaction among noise, rotor dynamics and 
characteristics of the optimization function as it applies to the control sys- 
tem stability and the total system performance. Consideration should also be 
given to implementing, the optimization function in ways that minimize its 
effect on control system stability. Such studies will also indicate the bene- 
fits that can' be obtained from signal conditioning and other forms of compen- 
sation for rotor and instrumentation effects that -would otherwise degrade the 
performance of the control system. These studies will involve the use of simu- 
lation models such as those developed for the feasibility study. Their objec- 
tive will be to develop the specific design criteria required for the develop- 
ment of a detail design. 

6.4 Rotor Characteristics 

The analysis of the control system. performance requires a method of repre- 
senting the characteristics of the rotor. This was achieved in the feasibility 
study by the surface coefficients of Table 13. The derivatives based on these 
coefficients should be verified to indicate that they represent the rotor to 
be used in the test configuration. Changes of these derivatives should be 
determined for the rotor with changes in rotor test conditions. This can most 
realistically be accomplished by comparing the characteristics as derived from 
the coefficients with characteristics derived from wind tunnel test data. 

Some test data have already been collected and are currently being analyzed, 
and additional wind tunnel testing is anticipated. The plans for this addi- 
tional testing should also consider the need to develop data appropriate to 
the analysis of the operation of the multicyclic control system. 
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DIMEN5IGM 21 200,16), X( 56) , Y( 56 J ,0 ( 56 ,56) ,p( 56 ,56} , 

2GI53 I, BI5bJ ,'ANS 1( 56) 

3, ,H'£ADl 20) ,r ITLESI 15, 5 ) 

‘CCM;<Gi\ /HDDAT'A/ KMOOEL 
COMMOK/T^O/N 

C0M.MJN/H0DELS/F( 200, 56), ZORIGI 200,16) 

5 FORMAT (IHl,//// , lOX, *SURGEN*=J==f=^ GENERATES A RESPONSE MATRIX TO FIT 
I A 5 ET OF DATA' POINTS*//) 


REAJ(1,505) HEAD 
50 5 F'jRMAT.IZD A4) 

C 

C 

C NDATPT = NO OF DATA-PGINTS (CASES) 
C NO OF COLJM'IS TO BE READ 


C 

C ■ ‘ ' 

C NMOD = 1' 

C ■ MODELl = M'JLTICYCLIC CONTROL THRU THIRD HARMONIC 
C • . NMOD =2 

C M0DEL2 = CTR DO, D1S,D1C 

C NMOD = 3 

C MODEL 3 - MULT ICYCL IC 4 HARMONICS 

e ■ 

R£A0( 1,13 ) NDATPT, NCOLS,NMOO, 1D8UG 
-NCTL=7 

IF( NMOD .EQ« 2) NCTL=3 
IF(NM00,EQ.3 ) NCTL = 6 . 

. 'IF{NM00.EQ.4i NCTL=9 '• 

ILEFT = NCOLS-NCTL 

R£A0( 1,5-10'H (TITLES! I,J ) ,J = 1, 5), 1=1, I LEFT) 

513-'F’oRMAT( 8-( 5A2)i . 

NVAR=43 

IF(NM00 ,EQ. 2JNVAR-14 
IF(NMOD ,EQ. 31NVAR=2,7 
TFINMOD, EQ.4)NV AR = 30 
IJ -FORMATIBIIO} 

INDEX=1 
K0UNT=0 

DO 20 1=1, NDATPT 
IFINMOD .NE.3J GO TO 22 

FUR MODEL 3 WANT TO SKIP FIRST THREE WORD ON CARD 00,D1S,D1C 
REA0(1,33) ( Z{ 1, J ),J=l,.NCOLSJ 
GO TO 20 

22 READ(1,33) ( Z ( I , J ) , J = 1, .N COL S ) 

23 - CONT INJE 
33 FORMAT (8F10. 2) 

33 FORMAT(33X,5Fi0.2/8Fl0.2) 

DO 25 1=1, NDATPT 
DO 25 J=l,NCOLS 
25 ZORIGI I,J J=Z I I, J ) 

Zl = NDATPr 
Z2<=NCOLS 



i'<i = ■< .' A .-ii- 1 
N^ = .NV AP.^-2 

35 50 re (5531t 5002, 5002,5004], NMOD 

5J0i CALL G?5rs {’ LOAD* , 'SURGMGOl* ) 

30 rO 5305 

5J0-2 CALL OPS irS(* LOAD* , 'SURGMOOi* J 
30 rO 5035 

5003 CALL OPSYS (* LOAD* , » SURGM0D3* } 

GO TO 5D05 

500 4. CALL OPSYS {* LOAD* SURGMOD4* ) 

500p continue 

00 31 I=1,NDATPT 
OQ. 31 J=1,N2 
51 F(I, J)=1.0 
W Rl T £ ( 3 » 5 ) 
ftRirE(3,5L6) 

51a format t//.) 

riRIT£{3-,505) HEAD 
GO TO (40,45,42, 48], NMOO 
4 3 CALL MODEL 1( ND'aTPT, INDEX] 

GO TO 46 • 

42 CALL MODEL3( NOATPT, INDEX ] 

GO TO 46 ' 

43 CALL MODEL2( NOATPT, INDEX ]■ 

GO TO 46 

43 CALL MOOEL4I NOATPT, INDEX • 

*5 KOUNT=KOJNT-H ' 

IF(NM00.EQ.2.AND.N.EQ-.8» GO TO 412 
DO 32 1=1, N1 
DO 32 J=1 ,N1 
Q(l, JJ=0«0 

32 P(I,J)=0.0 

DO 34 1=1, N1 , . • 

Y'{ I.]=0 iO 
-34 6U)=0.0 

DO 36 1=1, N2 
35 X(I)=1.0- 

DO 55 1=1, NOATPT 
DO 55 J=l,NCOLS 
55 Z( I , J] = ZORIG( I, J ) 

55=3.0 ■ 

S6=0.0 
S7 = 0 .0 

00 63 L=l, NOATPT 
DO 70 - 1=2, N2 
7 0 XU) = F(L,IJ 
00 83- 1=1, N1 
DO 93 J=1,N1 

93 P(I,J) = P( I,J)+X(I)#X( J ) 

33 Ym=Y(n+X( I)=<=X(N2} 

57=57+X(N2J 

53 S5 = 55+X( N2)=f=X( N2] 

Ml = P 1 1 , 1 } 

32=Y( 1 J 
DO 130 1=2, N1 
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J 



il=l -i 

si=pu, n 

53 = / ( I i 
S't = Pi 1,1) 

A=MI*S3-S 

5 8= ( MI *5^-Sl*Sl )4=( M 1=^=5 5-5 2^5 2} 

RAT I 0=0.3 • 

IFUA .EQ. 0.0) .AND. (BB .EQ. O.OJ) GO TO 61 
RATIO=A/SQRT(BB) 

51 R0UNO=nRATI0*lO00.+.5)/ 1000. ) 
rfRITE[3,110) II, ROUND 

UJ FCRMAH5X, » FACTOR* , 13, 5X, • CGRR ELATI0N = * , F6.3 ) 

U3--G( II )= A 

CALL OPSYSC *LQAD* ,»SURGMATX» J 
C.ALL HATINV'{P,Q, NG,Ni) 

IF( NG ,NE. 9). GO TO 105 
rtRir E(5,101) . 

IH F0R.MATI5X,10(«X:« J, • error IN MATRIX INVER SE *',1 0 (■=«=*) ) 

GO TO 412 ■ 

135 IF ( IDBJ'G .NE* 0 ) WRITE! 3, 102 ) ( { P ( NN ,MKr,NN=l ,N1 ) , MM=1 , N1 J 

IF { IDBJG- .NE. 0)WRITE( 3, 102JUQINN,MM),NN=1 ,N1) ,MM=1 ,NL) 

132 FORMAT!* MATRIX INVERSE * /6( 10F12 .4) // ) 

CALL MTRXHPI Nl, Nl, 1,Q,Y, Bt 0) 

WRITE(3,120) 

l>3 -FORMAT!/, 20X, ‘COEFFICIENTS OF BEST LINEAR FIT*/) 

DO 130 1=2, N1 

133 ANS11I)= U B! I )*1000. + .5J/1000.) 

WRl‘r!E(3,l40) (ANSI! IJ, I = 2,N1) 

14DF0RMAT(4!9F14.4/) ) 

C0NST1= ! ( B! IJ^IOOO.+ .SI/IOOO. } 

1.43 WRITE{3,150J CONSTl 

153 FORMAT ( 13X ,* CONST ANT TERH= * , F 10. 4// } “ 

WRirE!3,163) IT ITLESIKOUNT.J ),J = 1,5) 
lS3 FORMAT ( 23X,' TABLE OF RES lOUALS • / , 5X, ‘POI NT N0*-,7X,5A2) 

GOTO !18l,170,193, 194),NM0D 
17-3 WRITE(3,175J 

175 FORMAT !1H+,37X, • EST VALUE*, 4X, * DIFFERENC E • , 4X , * PE RCENT* ,4X,*DD* , 
14X,* DIS* , 3X, *D1C* / ) 

GO TO 187 

131 WRirE13,185) . ' • . 

GO TO 187 

135 FORMAT! 1H+,37X, * 6ST VALUE*, 4X, * DI FFERENS E 4X , * PE RCE NT* ,4X,*D0* , 
14X,* DIS* ,3X, * QIC* , 3X, * D2S* , 3X, * D2C* , 3X, *D3S * , 3X, *D3C* 

2,3X,*CQ * ,3X,*CLR/S'/ ) 

193 WRITEO, 1-92) , 

192 FORMAT !1H+, 37X, » EST VALUE* , 4X, * DIFFERENCE *, 4X , 'PERCENT* ,4X, 

1 ' »D2S*,3X, *02C*,3X, »D3S*,3X,*03C* 

2 ,5X, ’ D4S* , 3X,» 04C*/ J 
GO TO 187 
19'+ WRirE!3,198) 

193 FORMAT! 1H+,37X,*EST VALUE*, 4X, 'DIFFERENCE • ,4X, * PERCENT* ,4X,*D0* , 
14X,' DIS* , 3X, ‘Die, 3X, *D2S*-,3X, *D2C*,3X, *D3S*,3X, *D3C* 

2,3X,»D45* ,3X,*04C'/) 

137 12=3.0 

DO 130 J=1,NDATPT 
53=0.0 
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OQ 193 1=1, N1 
193 S8 = 58 + F( J, n*6( I) 

T2 = r 2+{F( J,N2)-S7/Zlj*=!=2 

S3=58-F(J,N23 

Z3=3 .0 

IF1F(J,N2J .EQ. 0.03 GO TO 191 
Z3 = S0^100./ F( J, N2 J 

I9i. ,^Rir E( 3,200 ) J, F{ J,N2), S8, S0,20, ( ZDR IG(J , JJ) , JJ = 1 ,NCTL) 
■lF(NM00.EQ,2 ) WRITE(3, 20 5 )ZGR IG( J , 8} , ZQR IG { J , 9) 

205 FORMAT (»+• ,90X, 2F8. 2) 

133 Z(J,1J=S0^ 

233 FORMAT H13 , F20. 3, 2F15.3, F 12.3, 5X, 9F 6. 2/) 

213 FORHAr(/////,20X, ‘ANALYSIS OF RESIDUALS*! 

NZERO= NOAT PT-1 
00 253- I=i,N2ER0 
H= I 

MM= I + 1 

DO 263- J=MM, NDATPT 
• IF(ZI M, U-ZU,i| J 260,260,270 
27 3 - M=J 
25 3 -CONT INJE 

IFIM .EQ. IJ GO TO 250 - 
PT=ZIM,1) , 

Z(M,1) = Z( 1,1) 
zu, i) = pr 
25 3 CONTINUE 
S=3^.3 
S2=3.3 

OU 280 -1=1, NDATPT 
S=S4-ZI 1,1 ) 

283 S2=S2 + Zn , 1)##2 
HI=5/Z1 

V=IZ1*S2-S=«‘S )/(Zl*l21-l. J J 
0=SQRT(V) 

WRirEt3,300) S,S2,MI,V,D 

303 FORMAT ( 23X,‘ SUM=‘ F15.4/9X, ‘SUM OF SQUARES=*P15.4/13X, ‘MEAN VALUE 
I F15.4,/,15X,*VARIANCE=*,F15.4/,5X,»STANDARD DEVI ATI 0N=* FI 5. 4/) 
DU = SQRT(IMI/B( IJ )=!==t2) 

IF(DJ-.OU 350,360,360 
353--Bm = B(l)-MI 

CONST 1= CONST 1-MI 
GO TO 145 

353 IF'U1.-S2/T2) ,LE. 0,0) GO TO 371 
R=5QRT( 1.-S2/T2) 

373 FORMAT (5X,« MULT CORR CO E FF=* F6 . 3/ ) 

WRir£(3,370) R 
371 WRITE<3,380). 

DO 403 1= 1, NDATPT 
RRR= I 

ZENO=(RRR-.5)=^100./Z1 
433 rtRITE<3,410) Z( 1,1), .ZENO 
41.3 = FORM AT { 2F15.3) 

412 IFIXOJNT-ILEFT) 35,420, 420 

333 FORMAT 1 23 X,' CUMULATIVE 0 ISTR IBUT ION * /lOX, ‘X * , 1 OX , *C UH PERCENT* 

;u /21X,' OF POPULATION'/) 

423 CALL EXIT 
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£N0 


/ ^ 

PHASE SJRGM0D1,=(= 

// EaEC EFORTRAN 

SJBRO'JTINE HQDELK NDATPT, INDEXI 
COMMON VMODAT A/ K 

COMMON/ M0QELS/f( 200t 56)»Z{ 200, 16) 



GO 

roti,2i, INDEX • 


'i 

K=3 


- 

2 

DO 

13 I=1,NDATPT 



FI I 

,2) = Zt 1,1) 



FM 

,3 J=Zt 1,2) 



FI I 

,4)=Zt 1,3) 



FI I 

,5)=ZtI,4) 



Fi 1 

,6 ) = Zt I, 5 ] 



FI I 

,7)=ZtI,6) ■ 


- 

Ft I 

,3)=Zt 1,7) 



F(I 

,9i=Z( I,1)-=S'*2 

• 


Ft 1 

, 10 ) = Zt I, iJ^Zt I, 

2) 


Ft I 

,ll) = Zt If l)*Zt I, 

3) 


Ft I,12J = ZU,l)*Zi I 

,4) 


Ft I 

, 13)=Zt 1, IJ’i'Zt I, 

5) 


Ft I 

,14) = Zt I, l)*Zt I, 

6) 


Ft I 

,15 )=Zt I, l)*Zt-I, 

7) 


Ft I 

, 16 ) = Zt I, 2)-^2 



Ft I 

♦ 17 ) = Zt I,2)^Zt I, 

3) 


Ft 1 

,18J=Zt I,2J*Zl It 

4) 


Ft I 

,19 ) = Zt I,2)*Zt I, 

5) 


Ft I 

,20 i = Zt i, 2)#Zt I, 

6) 


Ft I 

,21 )=Zt I,2)*Zt I, 

7) 


Ft I 

,22) = Zt I,.3)**2.- 



F t I 

,23)=Zt I,3)*Zt I, 

4) 


Ft I 

,24) = Z( i,3)?^Zt I, 

5) 


Ft I 

,25 )=Zt I,3)*Zt I, 

6) 


Ft I 

,26) = Zt I,3)#Zt I, 

7) 


Ft I 

,27) = Zt T,4)**2 



Ft I 

,23) = Zt I,4)*Zt I, 

5) 


Ft I 

,29) = Z( I,4)*Zt I, 

6) 


Ft I 

,30)=Zt I,4)=FZ( I, 

7) 


Ft I 

,31 ) = Z( I, 5)*=«'2 



Ft I 

,32)=Zt I, 5)*Z( I, 

6) 


Ft I 

,33) = Zt I, 5)*Zt I, 

7) 


Ft I 

, 34)=Z( I,6)=fr«2 



Ft I 

,35)=Z( I,6)^Zt I, 

7) 


Ft! 

,36) = ZtI,7)*=t=2 



.Ft I 

,37)=Z( I, 15)/ 10000. 


Ft I 

,33) = Z( I,15J^Z( I 

, D/IOOOO 


Ft I 

,39 )=Zt I, 15)«Zt I 

,2 )/lOODO 


Ft I 

,40 )=Zt I, 15)^Zt I 

,3)/l0000 


Ft I 

,41) = Zt I, 15)=f=Z( I 

, 4)/lCOOO 


Ft I 

,42)=Zt I, l5)*Zt I 

, 5)/10000 


Ft I 

,43)=Z( I, 15]^Z{ I 

, 6)/lCOOO 

LO 

Ft I 

,44) = Zt I, 15)=«=Z1 I 

, 7>/10000 


IN0ex=2 
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C it. - . iS THE NO. OF CONTROLS OF SERVE FLAP 

N “ 7 

JO 20 I=1,NDATPT 
20 F( I ,^5 ) = 7{ I,N1 
RE URN' 

END 

/ + 

phase S JRGMArXtSURGMODl 
// EaEC FFORTRAN 

SJBRQJTINE MATINV (A,B, . NG»NRA) 

, DIMENSION A(56,56), IR0W156), IC0L(56)»8( 56,56) 
NG=D 
N=NRA 
H= N+ 1 

DO 7 1=1, N 
lRGwm=I 
ICOLI n=i 
DO 7 J=1,N 

7 .B{ I, JJ = A( I, J ) 

DO 23 - K=l, N 
AMAX=B(K,K) 

DO 10 I=K,N 
DO 10 J=K,N 

IFIABSIBi I, J3 )-ABSC AMAX) )10,‘9,9 
9 AMAX=B{I,JJ 
IC=I 
JC=J 

U CONF INUE ■ 

KI«ICOL{IU 
ICOL( KJ=ICOL{ IC) 
iCOU ICJ=KI 
R I = I ROW I K ) 
i.RO‘W I Ki=IROW{ JCJ 
IROW { JC)=KI 
• IFIAMAXm, 12, 11 
12 NG=9 

GO TO Bo- 
ll DO 14 J=i,N 
E=B(K,JJ 
BIK, JJ = B( IC, J) 

14 B(IC,J)=E 
DO 15 1=1, N 
E=-B( I ,K) 

81 r,K) = BII,JC) 

15 B(I,JC)=E 
DO 16 1=1, N 
IFII-K)18,17, 18 

17 FUJ = 1 
GO ro 16 

18 F( I)=D- 
Id CONI INUE 

PVT= B(K,K) 

DO 3 J=1,N 

8 B( K, J) = B(R,J)/PVT 
F(K1 = F(KJ/PVT 

DO 19 1=1, N 


,F(56> 



IF{ I-K)2i,19,21 

21 AMULT=B{I,Ki 
00 22 J=1,N 

22 Bt I , J}=B( I, J ]-AHULT^B(K, J ) 

F( n = F-( n-AMULT=!=F(K J 

19 CONFINUS 

DO 20 1=1, N 

20 b{I,K) = F(U 
DO 25 I=i,N 
DO 24- L=1,N 

IF( IROW( D-U24, 23, 24 
2^ CONFINUE 

23 DO 25 J=1,N 
25 A( L, Ji = 3( I, J ) 

DO 2 5 J=’l,N 
JO 28 L=1,N 
iFiICGL{JJ-L)28, 29, 28 
2d 'CONIINUE 

29 DO 26 1=1, N 
25 Bn,.Lj = A( I,J ) 

30 RETJRN 
END 

SUBROJT INE MTRXMP(NRA,NCA,NCB, A, B,C ,ND lA G) 
DIMENSION AC 56,56), B( 56, 1J,CI 56, 1) 
IF(NDIAG)100,120, 140 
100 • DO 110 1=1, NRA 
DO irO J=1,NCA 
113 All , J)=-A( I, J) 

REIJRN 

120 DO 130 1=1, NRA 
DO 130 J=1,NCB 
CCI,J) = 0 
00 130 K=1,NCA 

13 0--CU ,-J) = C( I,J)+AC I,K)4=BCK,J } 

. RETJRN 

140 =D0 150 1=1, NRA 
DO 150 ^J = 1,NC8 
150 CCI, J) = A{I,I)=i=8C I,J) 

RETJRN 

-END 

PHASE SJRGH002,SURGM001 
EXEC FFQRTRAN 

SJBROUTINE M00EL2C NDATPT, INDEX) 

COMMON/ MODELS/ F( 200, 56), Z( 200, 16) 

COMMON /MODATA/ K 
COMMON/TWO/ N 
GO, T0(l,'2), INDEX 

1 K=0 = 

2 DO 10= 1=1, NDATPT 
FU,2)=Z(I,1) 

FCI,3J = ZU,2) 

F II ,4)=ZC 1,3 ) 

FCI,5)=Z(I,3) 

F( I,6) = 2C I,U*-2 
F( I,7J=ZC 1,2 
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FiiidJ = 2.iij3,---2 
Fli»9)=2( 

F( 1,10 ) = Z( I,l)'-Z{ 1,2) 

F( I., 11 ) = Z( I, i)^Z{ 1,3)’ 

Ff I, 12)=Z( I, 2)=«=Z( 1,3) 

Fi I,14) = Z{ l,9)‘ 

F( I ,15 )=Z( I,9)=J==f=2 
10 FU,13 J = Z( I, I, 2)=t=Z( 1,3) 

. IN0EX=2 
,(<- K+ 1 

C THREE IS THE NO, OF CONTROLS OF SERVO FLAP 

N=K+3 

IF1N.EQ,8) RCTJRN 
DO 20 I=1,NDATPT 
20 F{ 1,16 ) = Z( I,N1 . 

RETJRN 
• END 

/* 

PHASE SJRGM003,SURGMOD1 

// EXEC FFORTRAN- 

SJ'BROUTINE H0DEL3(’ NDATPT, INDEX) 

COMMON/ MOOELS/F( 200, 56), Z( 200, 16) 

COMMON /MODATA/ K »’ 

GO T0{ 1,'2) , INDEX 
1 K=3^ • . 

' 2 DO to 1=1, NDATPT 
‘F( I,2)=Z( 1,1) 

F1I,3)=ZI 1,2) 

FI 1 ,4) = Z{ 1,3) 

F( I,5) = Z(T,4) 

F{l,o) = Z{ 1,5) 

FI I,7) = Z( 1,6) 

Ft 1,8 )=Z( I, l)#-2 
F{I,9)==ZII,1)*ZU,2) 

F(I,10) = Zt I,l)*Zt 1,3) 

FtI,ll) = Z( I,li=^^Zt 1,4) 

Ft I,12) = Z( I,’i)^Z( I, 5) 

FtI,l3) = Zt I,1)*Z( 1,6) 

Ft I,14)=Z( I, 2)**2 
Ft I,15) = Zt I,2)*Zt 1,3) 

Ft I,16)=Zt 1, 2J^ZtI,4) 

Ft 1,17 )=Zt I, 2)=!'Zl I, 5) 

'Ft 1,1.8 )=Zt I, 2)*Zt I, 6) 

Ft I,19) = Z{ I,3J-=J^2 
Ft 1,20 ) = Z-{ I,3)4Z( 1,4) 

Ftl,21 )=Zt I, 3)*Zt I, 5) 

Ft I,22)=Z( IfSl’i'Zt 1,6) 

Ft I, 233 = Z( I, 4)=«'*2 
Ftl,24) = Zt I,4J=»=Zt I, 5) 

Ft I ,25 ) = Zt r, 4)=«=Z( I, 6) 

Ftl,26 ) = Zt I, 5)**2 
■ Ft I,27) = Zt I, 5)^Zt 1,6) 

13 -Ft 1,28 ) = Zt I, 6)*-2 
I.ND£X=2 
X= K+ 1 

C ' SIX IS NJM3ER OF SERVO FLAP CONTROLS 
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IM = K+6 

00 20 I=1,NDATPT 
23 F(i,291 = li l,N) 

RcTJRN 

END 

/ 

Phase 5 JRGM004,SURGMODl 
ff EXEC FFGRTRAN 

SUBROUT INE M0D6L4('NDATPT, INDEXJ 
COMMON/ MOOELS/F( 200, 56), Z( 200, 16) 
COMMON/ MO OAT A/K 
GO rO( 1,2), INDEX 
I K=0 


^ 00 

13 1=1, NDATPT 

F 1 1 

•2)=Z II, 1) 

F{ I 

,31 = 2.1 1,2) 

F( I 

,4)=ZII,3) 

F( I 

,51=21 1,4) 

FI I 

,6)=2I 1,5) 

FIX 

,7J = U 1,6) 

F( I 

,8) = 2(I,7)’ 

F(I 

,9 )=Z( I, 4)=i=*2 

FIX 

,13 )=2( I,4)*ZI 1,5) 

f tl 

, 11 )=ZII, 4)^Z( I, 6)- 

FIX 

,12) = Z( I,4)*ZI 1,7) 

FIX 

,13) = ZI I, 

FIX 

,14) = ZI I, 5)*Z( 1,6) 

FIX 

,15 )=ZI I, 5)*ZII,7) . 

FIX 

♦ i6) = Z( I, 6)44=2 

FIX 

,17 )=Z'( I, 6)4Zl 1,7) 

FI I 

, 18 )=ZI I, 7)442 

FIX 

,19 ) = ZI I, 8)442 

FI 1 

,20 ) = ZII,9).442 . 

FI I 

,21 ) = Z( I, 8)*Z( I, 9) 

FIX 

,22)=ZI I,4)4Zl I, 8) 

FIX 

,23)=Z( I, 4)4Z( X, 9) . 

FI I 

,24)=Zl I, 5J4ZI I, 8) 

FIX 

,25) = Z(I,5)4Z( I, 9) 

FIX 

,26 ) = 2I I,6)4ZI I, 8) 

FIX 

,27 )=ZI I, 6)4Z( I, 9) 

. FIX 

,23 } = ZI X,7J4Z( I, 8) 

F(I 

,29 )=ZI I,7)4Z( 1,9) 

FI I 

,33 ) = Z(I, 8) 

13 FIX 

,31 ) = ZI 1,9) 


INDEX = 2 
K.' = K.+ 1 

C NINE IS THE NUMBER OF SERVO FLAP CONTROLS 

N = K+9 • , 

00 20 I = 1, NDATPT 
23 FI 1,32) = 2( I,N J 
RETJRN 
END 
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AEROSPACE 

CORPORATION 


OLD WINDSOR ROAD. aUOOMFIELO. CONNECTICUT 06002 


KA\\A\ 

CORPORATION 

Report No. R-1494 
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APPENDIX B 
BASIC PROGRAM 


MCMODl 



KAAIAK CORPORATION 

OLD WINDSOR ROAD. BLOOMFIELD CONNECTICUT OSOOO 


KAittAK 

CORPORATION 


Mcnij: 

1 F.E 


e; if.oo 

r;u? Fil j fid' Fi3? " ^^‘5? R9 


Report No* R-1494 
21 January 1977 


£ij 

femb 

B2i? Bi ’ B£? B3? B4? E5^ B6^ B7? 18 

SO 

FERJ i 

r I;! j 1 • L ? C*3 ' C *4 ? C*5 ? L*t» ? » 



"III II ’i‘OU CHECI STflTEMEHT 

20 

IhPUT 

1 

2£ 

F i H 1 

"iriPUT THE LIMIT IMG nfiLUE 

65 

iliPUT 

11 

64 

PPIHT 

LIMiE.1 

65 

PPiriT 

“ii£S ii£C nss DSC ri4s H4C." 

62* 

FPIHT 



'r*0 

1^0 

110 

IcO 

i?0 

140 

1*50 

160 

170 

ISO 

c^lO 

£S0 

2S0 

£40 

£50 

£60 

£70 

£SO 

£S0 

300 

310 

3£0 

3£i 


330 

331 


3So 

331 

39£ 

333 

394 

395 

400 

401 
40£ 

403 

404 
500 
kii 
50£ 

503 

504 
600 
60 i 
60£ 

603 

604 
1000 


FOP :=-4 TO 4 
f]_r to 4 

FOR \ =-4 TO ‘i 
FOR L=-4 TO 4 
FOR r*1“-4 TO 4 
FOR \]=-4 TO 4 

31 “SOP r J) +30P I s +30R 

IF 31 : 3 THE^^ £40 

Z=fiO+Fi 1 I +fi£*^- J+fl3^t*.+fi4^L+R5-^r'1+fi6^r I I +fiS'> I I 

3=Z+E0^- 1 --^L+B 1 I 1 ^*H+E3^ J^-J+E4»‘ J^K>I!5-^' J^L+B6-^ J-^N 

Z=Z+B3^^K .j»L+CO^K^Tl+C 1 ^H+C £5>*L^L+C S-'^U^M+C 4^L-^H 

Z=Z+C 5^M=^M+C 
IF fiE3iZ:i *^0 THEM £40 
FFIWT USIMG £30U? 

INRGE 6{£D5£:tI^6H-£B 

r^E:4T M 

NEXT M 

NEXT L 

rCXT K 

r^EXT J 

ue:\t I 

PRIMT LIH(£) 

GOTO lEiOO 
PEN HP EQUrL 

BHTfi 744. 939? 1 . 3963? -1 . £54£^ -£ . 9£13^ -1 -'£5£S? - . £653? . 4£65? 1 . 5997? 1 . 1 £09 
BRTR * 173? - Fi714? 5684? • 175? 1 * 495? 0591 ?“- 019? 0E-£? 5E.94? *7733 

IiPlTfi .£434?-* 0533^ -.3934? -7101? “.0551? -3603? 1 . 3596? 0£95? 1. 1853 
REM ALPHA MAX EOUM. 

BATA 1 1 . 053 ^ . 1 96£ ^ . 5943 ? - - 0673 ^ - . 0093 ? - . 03 1 6 ? - . 1 30£ ? . 079£ ? - , 0£74 
BATA ".03£? .0503^ --03^“. 0149 1 . 1£49?“- 0613^ 0371 ^ .019^ -.0569? -041£ 

LATA — . UU J 1 « — - u5 56 ? • OiF,ci9 ? • 047 £ ? “• 086£? — . 0594 ? - o7u3? — . 00t»l ^ • U537 
REM BM EOUrL 

LATA 44 1 0 . 59 ^ - 1 4£ . ££ 3 ? -£06 - 873 ^ 3£- 5693? 5£. 894 ^ -1 0 . 43£7 


LATA 5 3 . 9363 ? 37 • 940 1 ' 3 . 6S7S ? —£0 - 5406 ^ — 1 £ . 9403 ? — £5 . 765 1 ? — ££ - 4££S 
BATA £9 . 9 1 4 1 ? - 1 £ . 5554 ^ -49 . 34£5 ^ -5 . 043£ ^ 35 . 07£3 ^ 60 . 6 1 77 ? 1 6 . 3733 
BATA - 1 3 - 5574 ? -93 - 5754 ? 1 6 - £7£7 ? 55 . 9666 ^ - 1 £ . 5085 ? 63 . 8309 ? 30- 8389 
BATA 37. 3327 
REM pi; EQUH. 

DATA £97 . 00 1 ? *- 1 7 . 5886 ? 9 . 3449 ^ £7 . 1 376 . 1 - 6 1 59 ? 48 . 9879 ? “ 1 5 . 7595 ? £ . 1 7 1 1 

BATA 8 - 8334 ? £ - 4898 ? -£ . 769 ? . 4 1 4 ? -7 . 7£3 ^ £ . 5843 ? £ - 5764 . 1 . 3003 ? -1 . 4£5 ^ - 1 4 . 4 1 04 


BAT A 6 . 075 ? 3 . 63 15? 16. 2^366 » — E* . 943£ ? 1 1 . 055 1 ? 1 *3 . 1 069 ' £4 - 1 043 ? 15. 1314?5. 659 1 
BATA 19.6047 
REM RH OP 

BAT A £37 . 3£ 1 ? — I £• . 5o93 ? 1 £ . 5836 ? £5 - 1 653? - 0t‘t«5 < 43 > 3534 ’—14- 46 1* 1 ' £ . 64US 
BATA 7 . 9£44 ? £ - 630 1 ? — £ • 7 £34 ^ 1 . 065 ? —7 .91 99 ? 1 £ . 32^08 ^ £ . 7869 » 2^ . 03 1 7 » — . 6995 
BATR “ 1 , y 1 9 1 ? 6 . 3857 ? 4 . £497 ^ 16*.! S03 ? —6 - 95E»6 ^11. 03 1 £ ? 10. 3£76 ? £ 1 . 4394 
BATA 1 4 . 758£? 6 . 1 £1 4 ? 1 9 . 1 995 
PEM RM IP 

BATA 79 . £904 ? -5 . 83£1 ^-6.1315-3. 8073 ? 3534 ? 3 . 7 1 7 1 ^ - 1 3. 5 1 48 ? . 535£ ^ i . 45£7 


BAT A . 0£37 ^ . 4007 ? — - 52*86 ^ 2* - O076 ? £ - 69J37 ? 1 • 2^302^ ' —6 . £94£ ^ —4 . 4792* ? — £ - 555 1 
BA ! A 1 . 3996 ? — - 4593 ? 1 . 1 62^6 ? — - 9 iJ97 ? 3 . 7933 ? 3 • 4.932» ^ 9 . U9 JF, ^ 10£9 

BATA -1.5£03?3. 1471 
Ei^B 
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APPENDIX C 


MAXIMUM SLOPE 


A flow diagram of the computer program used to investigate the maximum slope 
search characteristics is shown in Figure C-1 , and the listing is given in 
Table C-1. The first step is to set the initial values of the independent 
variables. These values are usually taken to correspond to no multicyclic con- 
trol, The increment of the independent variable to be used in the search is 
also set. 

The first step in the actual processing loop is to calculate the controlled 
variables (y's) using the coefficients in Table 11. Next, the value of the 
optimization parameter, P, and the derivatives of this function with respect 
to each of the controlled variables are calculated for the particular optimi- 
zation function being investigated. 

The decision on whether or not P has increased is bypassed during the first 
pass so that the next step is to calculate the derivatives of the controlled 
variables with respect to each independent variable using the coefficients of 
Table 13. The derivative of the optimization function with respect to each 
independent variable is then calculated as indicated in paragraph 3.1. If, at 
this point, the maximum derivative of the optimization function with respect 
to each of the independent variables is 0, the system has arrived at the "best" 
operating point and the process is stopped. If this condition is.not true, . 
the independent variable providing the maximum derivative of the optimization 
function is incremented and the process repeats from the step where the values 
of y were calculated. 

If the value of P does not decrease as a result of a particular iteration, it 
means that the operation has stepped across a valley of the surface describing 
the optimization parameter in terms of each of the independent variables. 
Therefore, the increment that caused that step is removed; the derivative of 
the optimization function with respect to that independent variable is set 
equal to 0; and the process is tried again using the independent variable with 
the next largest derivative. Note that if another independent variable causes 
a decrease in P, the whole process is repeated and a new derivative is calcu- 
lated for that which was previously set to 0. 

This process continues until an operating point is reached wherein increments 
of each independent variable causes the value of P to increase. At that point 
the process stops. However, the value of the increment can be changed and the 
operation continued from the operating point just found. 

The following outputs are printed for each iteration of the actual processing 
loop: 
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INITIALIZE x*s 
SET INCREMENT 


CALCULATE 

y‘s 


CALCULATE 
OPT PUNC, P 
AND EACH dP/dy 


DID P 
DECREASE? 


YES 

REMOVE INCREMENT 
SET CORRESPONDING 
SLOPE = 0 

1 



EA dy/dx AND 
J(dP/dy)(dy/dx) 


MAX SLOPE \ YES 
0 ? / 


INCREMENT x 
WITH MAX 
SLOPE 


Figure C-1* Shortest Path Model 
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i Ks table c-1. maximum slope model 


b 

13 

14 

15 
it. 

17 

■S' "S' 

■I'O 

CO 

£4 

£5 

£6 

3.1 

3£ 

oo 

oo 

3*5 


40 

41 
4£ 

43 

44 

46 

47 

48 
50 

t=:o 

._ic 

54 


IiflTFI 0 ” 0 > 0 ) 0 » 0 ■« 

MRT RERD 
FILES COEFF 
fll=-.5 
P£=1O00 
HIM ZC4]?YE4] 

IiflTfi 70005 750? 1£> 300 
MRT READ Z 
HIM RE6I 

HRTR 6"5 105 13s 15s 16s 16 
MRT. RERH R 
HIM MC4s3£3 
MRT RERH #isM 
IMRGE #5 311 
IMRGE #5X4H 
IMRGE #?4H 
GOTO 48 

PR I r IT " CHANGE "? 

•INPUT :=as:-?:£ 

IF X1=0 THEN 48 
IF S:i>6 THEN 47 
=>::■£ - 
GOTO 40 

INPUT XI 1 ] s :•=;[£] s :<i32 s j-:l 43 ? :=:c53 
REM 

FOR 1=1 TO 6 

PRINT USING 3i;a^<Lll 

ne:-:t I 


C63 


60 • GOSUE 500 
70 FOR 1=1 TO 4 

7£ PRINT USING 3£? V[ I 3+YCI f 1=3) 

74 nej=:t I 

76 GOSUB 6£0 

78 PRINT USING 3£s lO^sPl 

80 IF PKP£ THEN 90 

8£ :-=:EN3:=:4CN33.-fl£ 

84 SENS 3=0 

86 PRINT USING "#s£4X‘’ 

88 GOTO £5£ 

90 P£=P1 

100 FOR 1=1 TO 6 : 

iO£ SEI3=© 

110 FOR K=1 TO 4 

li£ Sl=MEKsI+13 

114 FOR J=1 TO 6 

1 i.6 R=M E K s I + J+R E I M I N J 3 3 

1 1-8 R=R®^=:E J3^ { 1+ C I=J) ) 

1£0 Si=Sl+fl 1 

1££ NEXT J , ! 

1£4 SEI3=SEI3+SiS^QEK3^-ZEK3 j 

130 NEXT ■ 

13£ PRINT USING 33s iO^SE 13 C-3 
140 NE:n I 

£5£ N3=T1=0 
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TABLE C-1, MAXIMUM SLOPE MODEL (continued) 

FDR 1-1 iO 6 ' 

if* FL'i'b I 3 .1 ii I HEi'1 fclScI 

258 Tl=flBS(SCI3) 

259 ,N3=I 

262 I 

270 IF H3=0 THEN 999 • 

290 n2=fl 1 ®SGH i S C H3 3,;i 

298 PRINT US IMG "2XSD"5N3S«SGm'R2J 

300 :'':CH33=:«;CH33-i-Fl2 

320 GOTO 50 

500 REM CRLC Y^S 

504 FOR K=i TO 4 

508 YCK3=MCKj13 

512 FOR 1=1 TO 6 

5 1 6 ‘f‘ C K 3 =‘r‘ L K 3 +M DK j 1 + n 5*:-^ C I ] 

52© FOR J=I TO 6 

524 Y C K 3 =’r‘ C F 3 +M E H U + J+R C 1 3 3 C 1 3 C J 3 

528 NEKT J 

532 he:==;t I 

544 he::=:t h 

548 RETURN 

820 REM CRLC PI S: DERIURTIUES ' 

821 IF W=4 THEN 681 

823 P 1 = 1 005K i; Y L 4 3 ^--Z C 4 3 J t2 
624 QC43=200:«Y[43.--ZT43 

629 IF W=3 THEN 885 
632 QC33=0 

834 -IF YC33<ZC33 THEN 639 

838 - P 1 =P 1 + 1 m i: 1 2* t 3 3 •••■Z E 3 3 - 1 J :i - T4 

638 Q C 3 3 =Q L3 3 +480® ( 1 2® E Y E 3 3 .-'Z C 3 3 - I D t3 ' 

839 IF W=2THEN 870 , - 

640 Pl=Pl+20®fYE23:''ZE23T-ir 

642 QC23=20 

644 IF YE23<ZC23 THEN 849 

648 Pl=Pl+l©®E30®i;YE23-^2E23-int4 

848 01 2 3 =0 C 2 3 + 1 200® f 30® C Y E 2 3 -^Z C 2 3 - 1 D t3 

649 IF W=1 THEN' 675 

850 Pl=Pl+36. 99+100® (YE 1 3-'-ZC 1 3-1 ;i 

852 QC 13=100 

854 IF Y[13<ZC13 THEN 660 

656 Pl=Pl + 10®i;7®i:YE13.'ZE13“l) )t4 

658 Q Cl 3 =Q E 1 3 +280® i‘ 7® ( Y C 1 3 .••-Z E 1 3 r 1 D t3 

680 • RETURN 

681 P1 = 1C10®YE43--ZE43 

662 C!E43=100 

865 P1=P1+100® (YE33-1 1 . 1 1 . ZE33 ’ 

866 - 0133=100 

670 . P'l =P 1 + 1 00® f Y E 2 3 -7453 .-’Z E 2 3 
67-1 0123=100 

8 75 . F 1 =P i F 1 00 ( ',V D -44 1 1 J • Z E 1.3 

876 >"ir 13=100 

680 - RLTURtI . . 

999 -.El III 
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Each independent variable 
Each controlled variable 
Optimization parameter 

Derivative of the optimization parameter with respect to each 
controlled variable 

Independent variable selected to be incremented and the direction. 
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A flow diagram of the computer program used to investigate the sequential 
search characteristics is shown in Figure D-1 , and the listing is given in 
Table D~l. The first step is to set the initial values of the independent 
variables and the slopes of the optimization function with respect to these 
variables. The independent variables are usually taken to correspond to no 
multicyclic pitch, and the initial slopes set to zero. The size of the incre- 
ment, the degree of lag and the amount of noise are also set as part of the , 
initialization process. 

The first step in the actual processing loop is to calculate the controlled 
variables (y's) using the coefficients in Table 11. Next, the effect of lag 
is introduced on the basis of the calculated values of the y's, the previous 
values of the y's, and the lag parameter set as part of the initialization pro- 
cess. These new values of y are stored for use in the next lag calculation. 
Next, the noise amplitude is generated for each of the y's and the values of y 
with lag and noise are used to calculate the optimization parameter, P, from 
the functions of Figure 3. It is the optimization parameter, affected by 
noise and lag, that is used in the subsequent control action. 

The decisions on whether or not P decreased and the basis for selecting the 
next independent variable to be- incremented are bypassed during the first pass 
so that the next step is to increment the first independent variable in the 
sequence. During subsequent passes, when P decreases, the change in P is used 
to calculate the slope of P to be associated with the most recently incre- 
mented independent variable. The basis for selecting the next independent 
variable to be incremented can then be alternated between a fixed sequence and 
the independent variable with the maximum slope. However, the maximum slope 
criterion is not used if the maximum slope is less than .01. 

If P did not decrease and if an increment of the opposite sign has already 
been tried, the increment is removed and a "restart" flag is set so that a new 
baseline value of P is determined prior to incrementing another independent 
variable. In addition, the sign of the slope associated with the most recently 
incremented independent variable is set equal to the smaller of the signs of 
the slopes determined with positive and negative increments. This aids the 
search process when the most recently incremented independent variable next 
comes up in this fixed search sequence. However, the magnitude of the slope 
is set less than .01, so that it does not become a candidate for incrementing 
in accordance with the maximum slope criterion. If the reverse increment has 
not been tried, and if the search strategy being investigated requires that it 
be done, the change in P is used to calculate and store a slope associated 
with the most recently incremented independent variable, the increment is 
removed, and one of opposite sign is applied. 
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Figure D-1 . Sequential Search Model. 
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7 MRT PEFiri 
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9 MfiT RERB S ' 

10 FILES COEFF 
11 . F1=0 

1£ F£=N1=H£=1 

1,3 Fli=-.5 

14 RESTORE 16 

15 DIM ZL4]jVi:4]jHC43 5EC43jDD43 

16 DRTR 7000 J 750 J 12 J 300 " 

17 MRT RERIi Z 
IS U— «00£5 

20 REr^rpRINT Fl: 0 ONE PI 5 1 NORMRLj £ ALL Y^S? 3 THREE PI’S. 

£1 REM F£: 0 RESTART, 1 1ST PASS’ £ -SEQ, 3 MR:-; SLOPE, 4 -INCREMENT 

40 P£=999 

41 DIM RE63 

4£ DRTR 6,10,13,15,16,16 
43- MRT READ R 

44 DIM MC4,3£3 

45 ‘ MRT READ #i?M 

46 DIM' R$Lll3,.B*-i:il3 

47 IMAGE #,3D 

48 IMAGE #,:-<4D 

49 IMAGE #,4D 

50 GOTO 58 

51 PRINT "CHANGE "3 

s£ input :-=:i,k£ 

53 IF. :-a*=0 THEM 58 

54 IF :«:i;-6 then 57 • 

55 :-:l:=:i 3 =:=<£ 


56 

57 

58 

59 

60 
61 
6 £ 

63 

64 

65 

66 

67 

68 
100 
101 


GUTU 
INPUT 
IF Fi=0 
PRINT ’■ 
PRINT “ 
PRINT " 
PRINT “ 


51 

Kl 13 , 1 -=Il £3 , i 


THEN 63 
£S £C 3S : 
OPT EM 
0PTN"5 

OC.O OCU o 


C 33 ,: 4 C 43 ,^=:i 53 ,:^^C '63 

3C 4S 4C"j 

HR 10RM RU"? 

3S S3C S4S S4C" 


C1=C£=C3= 


O' 


f- 4 =sC 5=:0 

IF Fl-0 THEN’ 100 
FOR I-l TO 6 
PR I, NT US I NG 47 j 1 3 
NEXT I 
REM 

IF F£>0 THEN 110 


103 MAT Y=B . 

105 GOTO 1£1 

110 .FOR K=1 TO'4" 

111 YCK3-MCK,U3 

li£ FOR 1=1 to- 6 ^ 

1 13. Y0K:3=:YtKiH^MLK, l.+ I*3^:^iC 1 3 
11-4 FOP J=I TO 6 ' ' 

1 i 5 Y‘ £K J -Y L K C ’ I +.J+R C 1 3 3 C 13 U 3 
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JL - 

li? 
120 
12i 
X dc! 

123 

124 

125 

126 
12 ? 
128 

130 

131 

132 

1 33 

134 
136 
138 
140 


J 

I 

r^E?;T ’ jabLE D-1. sequential search model (continued) 

Mfli I'-M-’ 

IF Fi^-0 THFr4 123 
GOSUd -!;gei 

.PRIf^T LiSIllG 48iPl 
IF FI <2 THEM 128 

PR I NT US I NG " 5D " 5 Y C 1 ] j 'f' C 2 3 f 1 [33? Y E 4 3 

PRINT USING "#j23?:" 

REM TIME LAG 
IF F2=i THEN 134 
MAT h=ci-t:i*h 
MAT Y=CTr=«‘!^ 

MAT Y=Y+H 
MAT H=Y 

IF F1=0 THEN 140 

PR I riT US I NG 48 ?Y C 1 3 5 Y C 2 3 j 1 G^Y E 3 3 j Y E 4 3 
GOSUE 600 


146 Cl=C:l+YE43 

147 C2=C2+YE43t2 

148 C3=C3+1 

149 C4=C4+P1 

150 C5=C5+Plt2 

165 REM GENERATES AND PRINTS NOISY DATA 

168 GOSUE 50E1 

169 IF FK2 THEN 180 

170 IF FI #3 THEN 172 

171 PRINT USING "#?4riK“jPl 

172 PRINT 


174 

176 


PRINT 

PRINT 


USING 

USING 


"#j 2-3fI" 

48 ! E 1 3 j Y E 2 3 ? 1 E 3 3 j Y E 4 3 


180 GOSUE 60© 

182 IF F 1=0 THEN 184 

183 PRINT " '■? 

184 PRINT USING "#s3D"jPi 

200 IF F2<2 THEN 221 

201 Sl=i:Pi-P2).-'A2 

202 IF PKP2 THEN 218 

203 PRINT " "P 

204 IF F2=4 THEN 210 

205 F2=4 

206 SEN33=S1 

207 N;CN33=>^EN33“2i«A2 

208 A2=-A2 

209 GOTO 300 

210 -^2=0 

212 SEN33=. 001»:SGN (RES (SEN33 ) 

213 :=-:CN33=:=aN33-R2 

214 IF F1=0' THEN 316 

215 PRINT . 

216 goto 316' 

218 SEN33=Sl 

219 PRINT 'ifl"; ■ 

220 GOTO 222 

221 PRiHT “E"? 

O I'O OO- 'C* 1 * 

U,L«L. r L,.“'r i 

223 Mi-! i' i—D, ■- 


MiH AESCSin 


p-4. 


KMIAK 

CORPORATION 

Report No. R-1494 
21 January 1977 



A 


OLD WINDSOR ROAD, BLOOMFIELD, CONNECTICUT ,06002 


Report No. R-V494 
21 January 1977 



CORPORATION 


c'^ij 


TABLE D-1. SEQUENTIAL SEARCH MODEL (continued) 

IF i ■■'HFH 


dDl 

£52 

£54 

£56 

£58 

£59 

£6£ 

£63' 

£64 

£70 

£71 


F£=3 

ri3=Tl=0 , 

FOR 1=1 TO 6 

IF fiBSOSniKTl THEN £6J 

Tl=fiBSCSm';i 

H3=l 

nej:t i 

IF TK.01 THEM £70 
GOTO 890 , 

N£=£ 

Ni=Sl+l-6s«CNl=6] 


C.I *i* 

N3=N1 



£90 

p£=flisi; 

SGN C S t H3 3 3 + C S C N3 3=0)1 

£91 

:=:cN33=:=: 

CN33+fl 

c. 

30G1 

IF F1=0 

THEN 

316 

304 

PRINT 

USING 

49? SC 1 3 s SC£3 5 SE33 ?SC43 r 3C5 

308 

PRINT 

USING 

”#j 311-'“ 5 F£ 

31£ 

PRINT 

USING 

":<ris'';H3»=sGNi:Fi£)- 

:31'3 

IF C3-;i 

8 THEM 

65 

3£0 

IF Fi#0 

THEM 

*“i F"i 

3£4 

PRINT 



3£5 

PRINT 

USING 

"3D. B" ? :«:c 1 3 5 :«:c£3-j y.z ~n ? nc43 

3£8 • 

PRINT 

USING 

" 4D " 5 C 1 --'03 ? SOR ( C£--'C3- ( C 1 ■- C 

•“'OO 

PRINT 

USING 

" 4D " 5 C4.-’C3.’ SQR (C5.-'C3- iC4^--C 

336 

GOTO 63 

1 


500 

REM flim 

i NOISE 


505 

FOR K=1 

1 0 4 


510 

Ui = l 



51£ 

U£=RND( 1.1 +.0008 

514 

IF U£<. 

500£ THEM 580 

516 

U£=U£-. 

5 


518 

Ul=-1 




5£0 

•5c!£ 

5£4 

'586 

530 

6£i0 

6£3 

634 

6*36 

640 

644 

t'46 

650 

654 


U3= 1 . 6£3 1 «SQR ( LOG (1 ."UE ) - . 3 J - 1 . 0 1 7 
U3=U3^«U1 

Y[K]=VtK3+U3SZnns«LI 

me:^:t k 

RETURN 

REM CfiLC OPT FUNCTION 
Pl=100ffii;YC4].-'ZC4 ] ) t£ 

IF VE33<2l3] THEM 640 

Pl=Pi+10»5 C !£:■« (YC3J.-ZC3]"1 1 ) t4 

P 1 =P 1 +£0^ ( Y C £ 3 . -Z L £ ] - 1 ;i 

IF YE£]:;ZE£3 THEN 650 

P 1 =P 1 + 1 0'^ i: 30^- C V i: £ 3 ••••Z C £ 3 - 1 ) J t4 

P 1 =P 1 +36 . 99+ 1 005K i; Y C 1 3 /-Z Cl 3 - 1 J 

IF YC13-;’ZC13 tkEM 660 

pi ■=rpi4 10ft t;7ft r vr 13 7 L 1 J“1 1 1 14 

*--ETURH 

ENIf 

D-5 





A 


iimm COR^PORATION 

OLD WINDSOR road. BLOOMFIELD. CONNECTICUT 06002 


KAAAX 

CORPORATION 


Report No. R-1494 
21 January 1977 


Once started, the model is allowed to run for a sufficient number of iterations 
to characterize its operation. Provisions are made for printing the following 
information on each iteration cycle: 

Each independent variable 

Each controlled variable as calculated, with lag, with lag and noise 

Optimization parameter from calculated controlled variables, with lag 
and with lag and noise 

The slope of the optimization parameter associated with each indepen- 
dent variable 

Independent variable selected to be incremented and the direction 

Mean and standard deviation about the mean for P and R^ over a number 
of iterations, usually 18. 
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